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ABSTRACT 


Rail track monitoring is very essential to have 
a smooth and safe ride on railcars* Standards have heen 
adopted to provide a minimum comfort level for a pas sen- 

• 

ger. Hence. it is necessary that a suitable track moni- 
toring system should be available for the measurement of 
track irregularities. 

At present midcliord measurement scheme has been 
adopted by Indian Railways* This has got poor frequency 
response in the pass band of interest* In this thesis, 
a better system, which malves use of accelerometers to 
pick up the track irregularities is proposed. Since the 
accelerometers give second derivatives of track irregula- 
rities, suitable digital filtering techniques using 
Digital Incremental Computation (DIG) scheme have been 
developed* 

Since track categorization will be done on the 
basis of power spectral measurement of the track irregu- 
larity, Fast Fourier Transform algorithms are used for 
the spectrum analysis. 

A completely independent dedicated signal pro- 
cessor consisting of an input digitizer, a control unit 
using 6800 microprocessor, memoiy modules for both pro- 
gram and data storage, cassette , RS 232 C console and 
teletype interface modules, a small console module con- 
sisting of fifteen digit eight segment LEDs and a Hex 



key board is developed 5 the system configuration and 
hardware details are described in detail. 

Tlie filtering and power spectral analysis is 
done on software basis on 6800,. Assembly language' "pro- 
grams are developed for this and are described in 
detail. 



CmPTER 1 


INTRODUCTION 

To have a smooth and safe ride of railway vehicles, 
the track tolerances should he within the acceptable limits 
and hence track monitoring is an essential part of railway 
operation* With the trend towards high speed transportation, 
the existing electrical and mechanical track irregularity 
measuring systems are inadequate and a better track monitor- 
ing system is proposed in this thesis. 

In this chapter various parameters of the track 
irregularities are described and the methods of measurement 
are briefly reviev/ed. Also an outline of the proposed system 
is presented* 

1.1 DEFINITION OF TRACK IRREGULARITIES 

Traditionally the track irregularities have been 
described by the following parameters* 

a) Unevenness of left rail 

b) Unevenness of right rail 

c) Alignment of left rail 

d) Alignment of right rail 

e) Gross level: 

f ) Gauge 

■ g)^ ; Twist y 1 • 
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However, the track geometry can be equivalently 
described by the following four parameters (see Fig, 1,l) 


a) 

Absolute centre 

line 

unevenness of the track - Z 

b) 

Absolute centre 

line 

alignment of the track - Y 

c) 

Cross level 


- XL 

d) 

Gauge 


- G 


Since the cross section of both rails change 
significantljr during the service, it is necessary to define 
these parameters unambiguously. A typical track geometry 
is shown in Fig.i,i, The various parameters are described 
below, 

a) Centre Line Unevenness of Track (z) 

The centre line unevenness of track is defined as 
the variation of midpoint of track in vertical plane 
from an arbitrary datum* This datum will be a horizontal 
plane for tangent track and standard changing horizontal 
plane (taking into account super elevations) in case of 
curved tracks (Fig, 1.,2) 

b) Centre line Alignment of Track (Y) 

The centre line alignment of track is defined as 
variation of midpoint of track in horizontal plane from 
an arbitrary datum. This datum is a straight line for 
tangent track and true circular or spiral shape in case 

of curved track (Fig. 1 * 3 ), 
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c) Cross Level (XL) 

If the tangent of two ra.il head crowns malce a 
slope V7ith horizontal, the cross level is defined as 
2L , 0 where 2L is normal spa.cing between the centre 
line of two rail heads,, and ^ is the angle the tangent 
makes with horizontal dattmi as showa in Fig. 1,1. 

d) Gauge (G) 

Gauge is measured as the distance between the two 
rail heads soa'ne distance 'h’ below the rail croim 
(h = l^-t- mm). 

e) Midpoint of the Track (M) 

The midpoint of the track is defined by first 
finding out midpoint of gauge line and then constructing 
a line perpendicular to the tangent of cro\ms of two 
rail head profiles. 

The traditional parameters can be derived from 
these basic parameters as follov/ss 

a) Unevenness of left rail will be Z + i^.L 

b) Unevenness of right rail will be Z - jZ),L 

c) Alignment of left rail will be Y + G/2 

d) Alignment of right rail will be Y - G/2 

e) Twist will be XL (1 ) - XL (2) / B where 

B = base, XL (l), XL (2) are the cross 


levels at the two end points of the base. 



In the track irregularity measurements, the loaded 
or unloaded conditions play a very important role and these 
irregularities should be measured only under loaded conditions. 
Considering the present day speeds, the deflections of 'the 
track due to its o\ra dynamics will be very small when compared 
to the loaded profile of the track geometry, hence are . 
ignored. Also the deflections due to track springing are also 
ignored. 

Since the moving railway vehicles will get the in- 
puts from the track geometry irregularities and depending 
upon the speed of the vehicle these inputs correspondingly 
excite the vehicle which results in vibration of the vehicle 
components. The railway vehicles are strong resonant systems 
hence this interaction is best understood in frequency domain. 
Each track geometry parameter is modelled as a random process 
and by finding out the direct and cross spectras of the ran- 
dom processes a more useful information can be obtained jji 
the design of track and Railway vehicles. 

1.2 EXISTING 1-IEASmEMENT SYSTEM 

Hie Indian Railways at present monitor the routes 
using the Amsler track recording cars which mainly monitor 
the track geometry parameters like vertical unevenness, 
gauge and twist. Such cars measure these parameters through 
’mid-chord versine measurement’ using a 3*6 metre chord. 

The length is chosen to be compatible with the existing 
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bogie length. From tlie track geometry parameters obtained 
through the ^*6 metre chord, the tracks are classified into 
the categories A, B, C, D, depending upon the tolerances. 
However, the wavelength of the track defects was completely 
ignored in this analysis as shown below. 

Let 'd' be the distance along the track and f(d) 
be the track irregularity. Let 'L< be the chord length used 
for the niidchord measurement scheme. Let AB be -.the drord 
and C is the midpoint of the chord. Considering the irregu- 
larity to be a pure sine wave form, different cases are 
illustrated in Fig, -1,4., From this it can be observed that 
the transfer function of the midchord at the ’peak inputs* 
can be approximated by 

2 2 

jn ( 5^)t = |(i - Cos I (1.1) 

where A is the wavelength of the irregularity. However, as 
the chord moves away from the peak (Fig. 1,4 case d), 
this approximation is not valid and the performance further 
deteriorates. Thus not only the value of X hut also the 
position of the chord affects the measurement of irregularity. 
Hence various frequencies present in the track irregularity 
are distorted by various amounts. 

For a dynamic system like a railway vehicle, the 
response is completely dependent upon the input frequency. 
Since the midchord measurement vitiates the track geometry 
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by attenuating the longer wavelengths (since for large A 
and for a given L, L « , the change in the irregularity 
is very small where as the irregularity with respect to a 
datum is not so ) , Amsler cars cannot measure the amplitudes 
of longer wavelengths. 

It is a necessary condition that the wavelengths 
■where the vehicles have substantial response should be least 
affected for the track irregularity measurements to have a 
relevance in estimating the vehicle response. Choosing such 
a chord which satisfies the above requirement for all vehicles 
at all speeds is a very difficult task. The speed variations 
come into picture because the vehicle frequencies in Hertz 
are related to the track irregularities through the speed of 
the vehicle, given by 

f = ^ (1.2) 

where 

V = speed of the vehicle in metres/sec. 

X = track irregularity wavelength in metres 
f = frequency in Hz, 

For higher speeds the larger wavelengths become 
more Important and the chord length for midchord measure- 
ments has to be longer to ensure that track irregularities 
for long, wavelengths are not attenuated. 
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A better system is proposed for tile track irregu- 
larity measurement which malces use of accelerometers mounted 
on axle boxes for picking up the irregularity waveform. Since 
the accelerometers have good response in the wavelength range 
of interest a better estimation of the irregularities is 
possible# 

1,3 PROPOSED SYSTEM 

In the proposed system accelerometers are placed 
on axle boxes (Fig,i,5^» Since the axle box is in contact 
with the rail v/heels, the accelerometers pick the random 
vibration of the wheel caused by the track irregularity. 

These signals are suitably processed to obtain the actual 
irregularity through the use of suitable signal processing 
methods. The power spectra of these signals are evaluated 
for the purpose of track categorization and in the design 
of tracks and railway vehicles. The measurement of rail 
vertical profile and lateral profile is given below. 

a) Absolute Vertical Profile 

To measure the absolute vertical profile, an 
accelerometer is mounted on each of the 2 axle boxes of one 
wheel axle set of the vehicle. Assuming that the wheels 
will always remain in contact with the rail in vertical 
plane, the output of the accelerometer will be double inte- 
grated and suitably filtered in the wavelength band of 2 m 
to ^0 m to acquire correct vertical profile. The central 
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line absolute -vertical profile of the track, the'ri, can be 
calculated by taking the average of the absolute vertical 
profile of left and right rail and the dynamic cross level 
will be evaluated by taking the difference of the absolute 
vertical profiles of left and right rails (Fig. 

b) Absolute Lateral Profile 

This measurement is also similar to the above. 

But because of -the standard play between the track gauge 
and wheel gauge, the assui-nption of wheel remaining in con- 
tact with the rail is no more valid in this case. Hence to 
measure absolute lateral profile, an accelerometer can be 
mounted on the axle box thus monitoring the lateral movement 
of the axle and distance of the inside rail face can be 
measured with respect to the axle with the help of two con- 
tact sensors, one each, looking into the rail face. Then 
the output of the two contact sensors and the lateral motions 
of the wheel set which will be obtained by double i ntegrating 
will be manipulated to obtain the absolute lateral profile 
of left and right rail (Fig, 1*6). Then the profiles are 
suitably filtered so that the correct information is availa- 
ble in the wavelength band of 2 m to 4-0 m and then by 
taking , the average of lateral profile of left and right 
rail will give -the central line alignment of "tie track and 
the difference gives the gauge of the track. Thus all the 
four basic parameters needed for the track geometry irre- 
gularity can be completely identified. 
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Data Processing 

After acquiring all the four track geometry 
parameters namely central line unevenness, central line 
alignment, cross level and gauge, are to he suitably proce- 
ssed for using the seune for track maintainence . For this 
processing, digital signal processing schemes will be much 
more useful, hence the data for every length of 512 points 
at, an interval of 2 m is proposed for analysing purpose. 

The power spectral density of each of the track parameter 
will be evaluated for this 512 points length and the value 
of roughness constant will be calculated by fitting the 
following model of Power Spectral Density (PSD) for the 
track irregularities. 

PSD = A i PSD is in mm^/cycle/M ; 

> is wavelength in M 
A is roughness constant. 

The value of A for each of the 4 track parameters can then 
be printed out, 

1.4 PROPOSED SIGNAL PROCESSOR 

As digital techniques offer various advantages 
over the analog schemes at these low frequencies, digital 
schemes are proposed- in this system. For this, the data 
from the accelerometer i s properly conditioned, digitized 
and then stored using appropriate signal conditioners and 
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Analog to Digital converters. The different operations 
are performed on the data to obtain the actual track irre- 
gularity. Although hardware digital filters can be used, 
as the power spectral measurement of the track irregulari- 
ties is also involved, instead of building different 
dedicated hardware supports the whole system is configured 
using a completely independent dedicated microcomputer with 
one of the easily available microprocessors and its family 
support. In this way, all the digital filter schemes, and 
power spectral measurement of the track irregularities can 
be obtained through the software developed for the micro- 
computer . 

6800 microprocessor and its family support is used 
in this system and a complete block diagram of the proposed 
system is shown in Fig, !+.l. The system consists of a 
central processing unit (CPU) module, memory modules for 
program and data storage, analog multiplexer and A/D con- 
verter modules for the sampling of input data and various 
interface modules to interact with the system. All these 
modules are assembled, tested and details are discussed in 
Chapter 4-. 
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OiAPTER 2 

FILTERING TECHNIQUES lOR TRACK IRREGULARITY TIEASUREIffiNT 


2.1 INTRODUCTION 

Tile proposed system for traxlc monitoring measures 
the track irregularities usi-ng accelerometers mounted on 
axle boxes of the rail carriage. The accelerometers pick 
up the second derivatives of the track irregultirities i.e., 

if a (d) is the irregularity, the output of the accelero- 

« « 

meter is a (d) where d is ttie distance along the track. 

There is thus a need for suitable filtei’ing techniques to 
obtain the actual track irregularity from the accelerometer 
output signal, 

Vdiether it is the measurement of vertical profile 
or lateral profile, an accelerometer ' is associated with each 
and therefore the method of obtaining the different irregu-* 
lari ties will be the same viz,, the output of the accelero- 
meter is to be suitably conditioned and then -double integra- 
ted, to obtain the actual irregularity a (d). Since the 
frequency of interest for the track irregularities is very lovr 
the design and stabili.ty of analog integrators at these fre- 
quencies pose various problems hence are not feasible for 
this application. A better approach will be to process by 
digital techniques in which the various inputs after samp- 
ling and storage are filtered using digital filters. These 
filters can be implemented either in hardvraro or software. 
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In the hardware scheme, the filters are implemented using 
arithmetic logic units, control and registers which are 
used for algebraic manipulation. In the software scheme, 
programs for various filters should be written which will 
then be processed using dedicated computers* 

In this chapter, digital filter schemes are des- 
cribed which can be used for this application. Assembly 
language programs for 6800 microprocessor applications 
have been developed and are listed in Appendix* 

2.2 DIGITAL INTEGRATORS 

As mentioned, to obtain the track irregularity, 
the accelerometer output signal should be double Integra- • 
grated. For adopting the digital schemes, these signals 
should be sampled and quantized using a proper Analog-to- 
Digital (A/D ) converter. Since the variation in the track 
irregularity cannot be very large with respect to distance, 
the frequencies of interest are very small. Hence a simple 
Simpson Integration for this application cad be adopted. 
This is sho^n in Fig, 2.1. 

Let f(x) be the accelerometer output; if f(x) 
is sampled at a constant spatial interval of ax, then 

k ^ ■’ :k, 

f (x) dx ^ ^ ^ ^^ 1 ^ 

^ ■ i=0 "■ '■ ■ ■ - ^ 

'' 

where Ax is the sampling distance and f (x^) is the 
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sampled and quantized value of the signal at distance 
Since .i& x, the sampling distance is constant, the inte- 
grator merely consists of an adder as given hy Eq, C2.l)j 
and a scaling factor. 

If this scheme is directly applied to the digital 
computer applications, Integrator may overflow due to the 
limited precision of the computers* Especially when a DC 
component is overriding on the signal, tliis DC level will 
be adding up causing the computer to overflow. Similarly, 
depending upon the input signal, the output of the 
Integrator may consist a DC_ level thus causing overflow 
for the second integrator to be followed which is used 
for the double integration purpose. Similarly very low 
frequencies also cause the overflow of integrators. 

Hence a high pass prefilter , which filters out all the DC 
and very low frequencies should be used before an integra- 
tor. A complete scheme for obtaining a track irregularity 
from the accelerometer is as shown in Pig, 2.2. In the 
following sections, the digital filtering schemes are 
described. 

2,3 DIGITAL FILTERS [1]. H 

Analogus to Analog filters, there are many 
types of digital filters like Butterworth, Tschebycheff , 
Elliptical etc. 1 ^ 23 . Hie digital filters can also be 
represented using zero-pole plots as in the case of Analog 
filters. Considering the real time application of track 
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profile measurement, the computation time for the ditigal 
filter should not be very high. And also the magnitude 
of the signal in the pass band of interest should not be 
affected by the filter response. Hence a single pole 
high pass Butterworth filter is chosen for this applica- 
tion although higher order filters give much sharper cut- 
off frequency. Since the processing time for Non-Recursive 
filters [23 is high, the recursive filter implementation 
is the best possible solution for this application. 

The transfer function of one-pole high-pass 
Butterworth filter in Z - domain is given by 


Y (Z) 

X (z7 


= H (Z) 



(w - 1 ) Z"** 


+ ¥ + 1 


(2.2) 


where 

Z = exp (2TTi f/N) 
w = Tan ( TT f /N ) 

N = Number of samples per second 
f = Cut off frequency of the filter in cycles/sec. 

f = any frequency for which transfer function is 

to be calculated. 
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From this equation, the recursive relation for the filter 
is 

h = (tIv) 0^ - ^i-1 h - h-i] '2-3> 

where f^^ represents the present value, is 

state earlier value. represents the present input and 
is the present output. The above recursive equation 
can be realized as given in Fig, (2.3)* 

As it can be seen from this realization, since 
the present input is subtracted from tine previous input, 
the DC component is removed at the output of the filter. 
Since w is a function of cut off frequency, by properly ■ 
choosing a value for w, the cut off frequency can be fixed* 
Since all the computer operations will be in binary arith- 
metic, it is better to choose a ’w’ which is in multiples 
of 2 so that divisions and multiplications can be performed 
very easily. In the practical track monitoring routines, 
it was observed that a cut off frequency of 0,12? Hz is 
suitable .Considering the present speed range of Interest 
(12O KMPH) and sampling rate of 5l2 samples per every Km, 
a value of 2 for w is chosen. 

For the integration, the trapezoidal integra- 
tion scheme also works out to be better, for which the 
transfer function is given by 


( 





I (Z) 


k 

2 


) 


(2.^) 
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from which the recursive relation works out to be 

X. + X,. , : 

^i = ^i-1 ^ ^ ^ ^2.5') 

where is the present value, f is one state earlier 
value, h is the sampling interval. The recursive relation 
for high pass filter and a trapezoid integrator combined 
together will be 

^i ~ ^ 1 + w ^ "^i-l ^ “ 2 ^ ^ ^ 

( 2 . 6 ) 

The composite filter is simulated on DEC 1090 
computer. Whether it is the Simpson or trapezoidal inte- 
gration , the following observations are made for the 
composite filter 

a) Settling time of the composite filter is high. 

b) \^^len the complete double integration scheme is 
implemented, due to the settling time error, 
tiiere is a fluctuating DC component observed at 
the output • 

As this setting time is large, the filtered 
output doesn’t correspond to the actual profile for this 
period. Hence, during this time the data is to be ignored. 
A better alternative in which settling time can be reduced 
is the ’incremental computation scheme ' which works on 
the principle of differential input - output scheme in 
which increments in the data are used for processing. 
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2.4 INCREMENTAL COMPUTATION SCHEME EOR 
DIGITAL FILTERING OF LOW FREQUENCY 
SIGNALS C43 

Normally, for the digital filters which operate 
at low frequencies, the ratio of cut off to sampling 
frequency is very small. In other words, for the filters 
which have low ratios of cut-off to sampling frequency, 
the poles are located close to the unit circle in the 
Z - plane. The recursive filters which operate with the 
poles very close to Z = 1 and use finite word lengths 
are prone to the effects of sensitivity and round off 
errors. A great deal of work has been done to reduce 
these errors and also different schemes have been presen- 
ted in using the incremental computation for the low 
frequency signals. Detailed discussions can be found in 
the reference [VJ. 

2.4.1 Difference Between Normal and 
Incremental Computation £4} 

Definitions 

a) Normal Computation s A computation procedure 
that evaluates each value the dependent variable 
assumes independent of its previous values is 
referred to as the normal computation scheme. 

The value of the dependent variable computed in 
this manner is termed ‘normal value ’ • 



25 


b) Incremental Computations Normally many physical 
system are slowly' varying. In ^ther words, for 
a given change in the independent variable, the 
amount by which the dependent variable changes 
from its previous value is very small* Bie 
succeeding value that the dependent variable 
assumes can be computed by adding an increment 
or decrement to the one previously computed. 

A scheme that employs this principle in evalua- 
ting values that the dependent variable assimes 
is referred to as an incremental computation. 

2*h-m2 Using Digital Incremental Computation (DIC) 
for Digital Filtering 

Let the required transfer function be 
H(Z)=:y(Z)/X(Z). In thu conventional filters where 
the normal scheme of compnbation is adoptr-d, the input 
and output are not in incremental form. In the DIC scheme 
the input is followed by a differencer, a transfer func- 
tion H' Cz) for which the input and output are in incre- 
mental form and an output integrator as shown in Fig. 

As from this illustration, for the two structures to be 
identical the transfer functions H (z) and H» (Z) should 
be identical. A simple DIG structure of first order filter 
is showi in Fig. 2.5« ©-j (n) and Cn) are the noise 

sources arising from the quantization operations for multi- 
plication and integration* The transfer function H’ (z) 
of the DIC structure is 
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H> (Z) 


Y (Z) 

X '(z)' 




YlI 


B Z 
o 


•1 


■1 2 "' 


( 2 . 7 ) 


The transfer function of the first order conventional 
filter is 


H (Z) = 


Y (Z) 

x“Tz7 


c - d Z 


-1 


cl Zj 


7-* i. 


( 2 . 8 ) 


In order to the transfer functions to be equivalent 
corresponding coefficients in transfer functions are to 
be equated. From this, the DIG filter coefficients are 
given by 


Bo = d 

=: (c - d) / (2.9) 

= (1 - a) / G^ 

where G^ is the gain of the integrator. 

A block diagram of the DIG filter with inputs 

and outputs in normal form is shown in Fig, 2.6. The 

coefficients A^, will be chosen such that E' Cz) = 

G H (Z) where G is a scale factor and is imposed 
o o 

by the structure since it operates upon the increments 
of the input. And also Go can be used for controlling the 
dynamic range of the filter. Normally Go should be chosen 

such that 4^'^! jnax “ max* ^ above DIG 

structured filter, the input and output are in incremental 
form. In order to specify the complete structure of the 
DIG, the scheme in Fig, 2.6 will be adopted where input 
and output are in normal form. 
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Second Order .Filter 

A second order filter can also be derived using 
DIG schemes on similar lines* Consider the second order 
transfer function 


H (Z) = (2.10) 

1 - a Z ' + b Z“^ 

The structures of second order conventional and DIG fil- 
ters are given in Fig. 2.7 and Fig, 2.8. From Fig. 2.8 
the transfer function of DIG filter is given by 
, (B - B G + B, G G ) - (2 B + B G ) B^ Z^^ 

g’ (2) _ 2 J 2 iS 1__2 2 J 2 

1 - (2 - G^ - A^ G^ G^) Z"^ + (1 - A^ G^) Z"^ 

( 2 . 11 ) 


The corresponding coefficients are 

B, = (d - 2) / 0„ 

B^ = (l + c - d) / G^ G^ 

h = <1 - 

— (l “ a + b) / G_^ G^ 

-b. 

where G^ = 2 is the gain of the integrator, 

e^ (n), (n), e^ (n) are the noise sources arising from 

the quantization operation of multiplication and integra- 
tion, 

2.5 NOISE STITDYj DESIGN ESOGEDURe M 

The roundoff errors that occur in recursive 
digital filters employing fixed point arithmetic are mainly 
due to the finite world length of the digital computers. 
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The first error is caused by the addition of two fixed- 
point numbers resulting i. overflow. By restricting the 
dynamic range of the signal level, the overflow can be 
controlled not only at the input and output but also at 
other intermediate points. The second error is due to the 
multiplication of the signal with the filter coefficients. 
When two numbers are multiplied, the resulting product 
word length is not equal to the sum of the lengths of tbe 
two words, being multiplied. In a recursive filter, this 
product is again multiplied vrith a sample of the input 
signal in the next recursion. But the length of the pro- 
duct should be quantized to limit the register length. 

This quantization of the product word gives rise to an 
error that accumulates as succeeding multiplications are 
made. The model for fixed point round off noise following 
a multiplication is shown in Fig. 2*9. 

2.5»1 Koise Power in Second Order 
Recursive Filters 

Consider a transfer function for second order 
recursive filter 

1 - a Z" ' + b 

For obtaining the model of recursive filter 
the following assumptions are made. The operation of 
quantization used is rounding. 
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a) Any two different samples from the same noise 
sources are uncorrelated. 

b) Any two noise sources associated with different 
multipliers are uncorrelated. 

c) Each noise source is uncorrelated with the input 
sequence. 

Let the variable be represented by (b ’ + 1 ) 
bits in the fixed point arithmetic which includes the 
sign bit also. Then each noise source can be modelled as 
a discrete stationary white random process with uniform 
power spectral density as shown in Eig. 2*10, 


With the above assumptions the filter can be 
modelled as shown in Fig. 2.11. When the input signal is 
completely absent or becomes constant, the three assump- 
tions of uncorr elatedness fail and the filter exhibits an 
oscillatory behaviour at tne output which is termed 
LIMIT CYCLE OSCILLATIOH. The round off noise power at 
the output is given by M 


i -f- 13 

_ g ' -.‘' 13 ) (-1 + b - a) (1 + b 4 - a) 


for the conventional filter. 


(2.1^) 



A similar expression can be derived for the DIC filter 
given by 

„ cf (1 + b) + {3 + 2 (3 “ a - t*! (1 - a + b) 

b, 2 -...-..-Y... iL.. 

__ !' 6 (1 «• b) vl - a + b) (1 + a + b) 


( 2 . 15 ) 


vjiiere =2 axid = 2 are the gains 

•-7 

of the Integrator. Normally is of the order of 2 

P 

hence can he neglected in the above equation. 


CTaS. 0:1 — ^ - a _ j:,h ^ 

3 (1 -h)(l + a +h) 


or 


®d2 ^ 


2 4 - G ^ 


3 A^ ;(4 -'g G^ G^; 


(2.16) 


From this it can be seen that G^ and G^ affect 
the round off errors. The expressions for and A2 are 
(from 2.12) 

lo 

A^ = (1 - b) / G^ = (1 - b) 2 2 ( 2 . 17 ) 

■b- 6 

A^ = (1 - a + b) / G^ G^ = (1 - a + b) 2 . 2 

(2.18) 

Choice of b^ • 

From the above it can be seen that A^ depends 
only on G^ while A^. depends upon both G^ and 62. There- 
fore depending upon the values assigned for G^ and G2 in 
the numerator of (2.I6), the follovong situations will 
arise 

a) 2 A|« A2 G^ ( 2 . 19 ) 

b ) 2 A^ » >2 ^ (2.2O) 

and in the denominator, normally, 
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Using (2.19), (2.21 ) In (2.16) 

/ 12 ( 2 . 22 ) 

From this, small values of A^ should be avoided, 
otherwise round off error increases. It was shown []^3» 
from experimental studies, that values of A^ < 1 should 
be avoided. Therefore A^ ^ 1 is preferred to minimize 
the round off error. Using this condition in (2.17), 

(1 - b) 2 ^ ^1 
or . 

^2 5“ ^°®2 ^ 1~^b ^ (2.23 ) 


and the value of G2 can be obtained by choosing a value 
for b^ using the above equation. 

Choice of b^ s 


A choice for b^ can be obtained by setting 

mif ^ in Fig, 2.8 . Let (n) and 

h^ (n) be the impulse responses from the input nodes to 
those at which y^ (n) and y^ (n) present 


Then 


|y,ji 


- y 2^ |h' Cn)j 




max 


max I ^2 


(2.24) 

(2.25) 


where 
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(Z) 


_ (1 + c - d) 4 - (ad -- a-d-bc + 1 ) z"'^ + (i-a-A*) Z*"^ 


rw 


V 

n 


h,' (n) 


-n 


(2.26) 


(Z) 


(c -- 1) +(a--d) 7 r\ + (1 --13) tT'^ 
G^ D (Z) 

= 5 (n) Z"^ 


n 


D (Z) 


1 - a Z”'* + b Z""^ 


(2.27) 

(2.28) 


Using these equations, 


^1 ^2 ^ ^2 n" i ^2 


or 


or 




1 I 


n 


jh^ ta)lc-| Ih^ (n)| 

^ |h-| fa)j 


|K(n)l 


(2.29) 


Choice for dynamic range scaling factor G^'. 

From the Fig. (2.6), v^rhere the implementation 
of DIG filter using the input and output are in normal 
form is presented, the scaling factor G^^ can be obtained 

max ~ ^ ^ max 


by choosing j-Ayj or from the figure 


it is the same as 


(1,1 


= t^t 


max 


max 


C2.3O) 


Using this G^ can be evaluated as 


G _ ^ _ 3 __Sax 


max f b <n)| 


(2.31) 
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where h (n) is the impulse response of the filter. 


In this ^ |h (n)| is the infinite summation and a 
n 

reasonable value of n should be chosen such that the 
error in the siimmation is small. In the extreme case 




max 





h (n)| 


(2.32) 


The value of obtained from this equation 

can be adjusted to the nearest binary number in order to 
eliminate the multiplier • FORTRAN programs for calcula~ 
tion of these parameters are given in fV} » 


2*6 HIGH PASS FILTERING USING DIG SaiEME 

The transfer function for a second order high 


pass filter is of the form 

2 

f r7~ I ' 

H CZ) 


(1 - z-U 

.-1 


1 “ a 


+ b Z 


-2 


C2.33) 


Comparing this with the transfer function of the second 
order filter 


H (Z) = 


d Z”'* + Z'2 


1 - a Z""^ + b Z"^ 


(2.3^) 


we get e=1? d=2; using this the coefficients for 

DIG structure can be obtained as 


B 


= 1 


= Op 


.Bp - ^ 


(2.35) 


= (1 - b) / G^P A^ = (1 - a +b) / G^ G 


1 2 


The block diagram is shown in Fig. 2 •12. 
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2.7 high pass filter for track profile 

MEASUREMiATT 

As described earlier, a one pole Butterwortli 
high pass filter is first implemented for the real time 
application of track profile measurement. For this the 
high pass function considered will he 

H CZ) = (2.36) 

(w-l)Z +W+1 

For 'Chis filter, the working programs in 4- , 
are im-plemented on DEC- 1 090 computer and the following 
coefficients for the DIG scheme are obtained. 

= 1 = 1 1.2 = 5 = 1 

and the DIC high pass filter will be as shown in Fig. 

2.13. 

The proposed system thus would consist of a 
first order Butterworth high pass filter using DIC 
structure and a simpson integrator. An assembly lan- 
guage program for 6800 microprocessor has been developed 
and is given in Appendix A. This program is for the 
composite filter consisting of a high pass filter and 
integrator. The programme is tested on 6800 micro- 
processor development system and the results and discu- 
ssions are presented in Chapter (5). 

After obtaining the different track irregula- 
rity parameters namely vertlcle profile of left and Tight 
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rails, throug’i simple a.lgebraic Eianipulation described 
in Chapter 1, one can derive all the other irregulari- 
ties like unevenness, ■ gauge etc. 

In th.e following Chapter, the power spectral 
measurement of these irregularities through Fast Fourier 
Transform (FFT) technique is described. 
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CHAPTER 3 


POVvTHl SPEGTRJIL MEASUREMENT OF TRACK IRREGULARITIES 
3*1 INTRODUCTION 

As tile trains move on irregular tracks the 
comfort of the passenger is affected. Standards have 
been established for minimum comfort to be provided. 
Equivalently standards are also laid down for categori- 
zation of rail tracks through roughness factors. These 
factors are intimately connected to the Power Spectral 
Density (PSD) of the random process correspondent to 
track irregularities. 

In Chapter 2 the methods of obtaining the 
waveforms corresponding to the track irregularities are 
discussed. Using this the power spectral density asso- 
ciated with the irregularity is to be obtained. A set 
of N samples associated \^ith the track irregularity 
is available to estimate the power spectral density. 

In the following sections, the estimation of PSD through 
the use of Fast Fourier Transform (FFT) is discussed. 

The implementation of FFT on 6800 microprocessor is also 
described. 

3.2 POWER SPECTRAL DENSITY (PSD) 

The random process is said to be *strongly 
stationary* \&ien all the possible laoments of the process 
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are time invariant. The process is said to be 'weakly 
stationarj^’ in vjide sense if the first two moments namely 
mean and autocorrelation of the process are time invariant. 
Practically it was observed that ferack irregularity 
parajneters constitute a weakly stationary random process. 

The PSD 'of a process is defined by 

1X3 

(<^) = J (7 ) e" dr (3,1) 

00 

where X is the random process, S^ is the PSD 

and R^ ( 7 ) is the autocorrelation of X, The PSD is 
non negative for all OJ and even function. Because of 
this the PSD is defined only over the positive frequen- 
cies from Zero to +00 and is referred to as 'one 
sided spectrum'. The ordinates of one sided spectrum 
will be double the right half of mathematically genera- 
ted two sided spectral density, and is used here. The 
important constraint underlying the definition of PSD 
is that the sample function (or realization) of the pro- 
cess should be a record from a stationary random process 
of infinity. However, the concept of PSD is also equally 
applicable to weaic stationary process. Since it is not 
possible to take a record length of infinity, a realiza- 
ble length L will be chosen to make the definition of 
PSD physically realizable and to obtain a consistent 
estimate of PSD. . 
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3*2.1 I-lathematical Analysis {"51 

Although the concept of randora process has been 
developed for continuous -time random signals, 41ie same 
can be easily applied to the case of discrete-time random 
signals* The mean and autocorrelation of a discrete 
process X (nT) can be expressed as 


03 


E j 

A (nT)3 

= Jx p (x 5 nT) dx = 

(3.2) 






(kT) 

= E [X (nT) X (nT + kT^ 

(3.3) 

or Rv 

A 

(k) 

= E fx X , 1 

L n n4c-l 

(3.4) 

where k 

is the 

time separation or lag. 



In the discrete case, the N samples from a 
realisation of a zero-mean stationary random process 
will be used to ’estimate' its autocorrelation function, 
which is defined as 


H (k) = E J (3.5) 

First the estimation R (0) = will be 
considered. Let the estimate be 

N-1 ^ 


ys 2 

CT 


* 2 2 2 

_ R (0) - A y 

_ ^ ^ x^ _ 


In this, if k is large x^ and x. 


14c 


(3.6) 

will be effec- 


tively independent eventhough the realizations x^ , 
o 

x| , ... need not be independent. 
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The estimate is eq. (3*6 ) is unhiased 

(for large N) since = pr^ and also consistent 

since lim E f (cr ^ - Cr^ >3=0 i.e. variance is 
N-> oo ^ 

zero. Hence for large N, the ^ from a single realiza- 
tion has a higti probability of nearing the true value 


In the similar lines, the estimate for Ck), 
the autocorrelation function at other lags can be pro- 
posed as 


Rjj (k) 


^o ^>kt ^ ^tki ^ 1 + 

N 


1 

N 


H -}ki. - 1 

V 

iC- X. 

i=0 ^ 


^i +!kj 


(3.7) 

k = Oj + "J j + 2,... 

(3.8) 


+ N-1 


v/hich averages all possible products of samples separated 
by a lag of k and for large k .there are very few possible 
products. For k N, the re are no possible pairs with 
this lag available, hence it will bo assumed that the 
estimate of Rt,t (k) is zero. 

iI-Wci-1 

^0 ^ = °> ■±’'> i “-I 

0 ■ k ^ H 

(3.9) 

(k) is an evonfunction of k, as is true in the case 
of R (k). Here N denotes that the^ realization is of 
N samples in length* 
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The mean Talue of the above estimate is 

r- — \ ^ N— 1 •• i h I 

E lR„ (k)J = u 2: ® =" 3 +)kt'] 

3=0 

= ( 1 - ) R (k) (3..IO) 

Thus the sample autocorrelation function in (3.9) is 
biased since the mean of the estimate is not the true 
autocorrelation function R (k) at lag k. However, as 
the term N vanishes. Hence Rj^ (k) is 

assymptotically unbiased. An unbiased estimate can be 
easily obtained by dividing the sum in (3.10) by (N - |kl ) 
rather than N, Then 

H-l_-fki 

® 1% = N - jkl itn ® ^ 3 * 11 ) 

Hov/ever it was shown [53 "that for the estimate 
of (3*10) variance is small for N » k lor which bias 
is also small. Therefore, unless a large number of lag 
products is available for an estimate at a k, the 
resulting estimate Rj^ (k) will not be reliable. 

3.2.2 Estimates of PSD 

The power spectral density and autocorrelation 

function of a realization are related through the 

Discrete Fourier Transform (DPT), given by 

00 

= 1 R (k) 

k=: - 00 , ■ 


s (ca) 


(3*12) 



k -7 


Considering the estimate of R (k) as (k), given 
by (3*9), the spectral estimate (00) ^ill he, 


Sjj (^) 


CO 

X Rjj (k) e’ 

K=: -CO 


(3.13) 


The spectral estimate S^j (c^) is also related to the 


Sn = f 


N“1 -0«-oT j 2 

X ^ e 
n =:0 ^ 


data samples (00) through the DFT, 

TVr ^ -1 *2 ^ m ^ 

'i = i 

(3*1^) 

This can he proved as follows. Consider the truncated 
signal 




0 


Then 


Rjj (k) 


1 

N 


n = 0 , 1 , ... N - 1 
otherwise 


00 

1 

n = -CO 


N N 


( 3 . 15 ) 


(3.16) 


From (3.13) Sjj (CO) is given hy 


00 
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T 


k= -CO n = -00 


N jj -jkcoT 
^ ^4k ® 


^ n= -co^ k= -00 


® (3.17) 

-jco (n 4 c)T 


= = i |Xh <“)| 


(3.18) 
2 

(3.19) 


CO 


in which Xjj ( 03 ) 


Z . 

n= -00 


jj -jncoTg N-1 


-jnOoT^ 

E ® 

n =0 " 
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Considering the estimate (GJ), the mean of ('^ ) is, 


N 


oo 


ELSj, (CO)] = _ r &)] 


rjkCA>T 


k=: -oo 


e 


(3.20) 


From (3,10), 

E (to)] 


N 

z 

k=-N 


R (k) ( 1 - ‘ 


IM 


^ -jk<A>T 
) e ■ ^ 

C3«21) 


liius the spectral estimate is the biased esti- 
mate of S f to), which differs from S by the term 




V, 


N 


(1 „ iliL) 


k 


0 


Ikt ^ N 
|ki > N 


(3*22) 


However the (Ob) is assymptotically unbiased since 
the term |kf / N vanishes for large R. The term v^ is 
called the ’window function ’ • The presence of this 
term in spectral estimate causes the power that presents 
in a particular frequency range to be smeared out over a 
wider range. This spreading of the power is also called 
as ’LEAKA.GE*. For having an accurate estimate of PSD, 
suppression of sidelobes or reduction of leakage is 
essential. This can be achieved by using spectral win- 
dows on raw power spectral density. There are many 
window functions available [53 3^^ ihe following are the 
most popular ^^^indows that are practically used. 

a) Rectangular 

b) Hann 

c ) Hamming 
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d) Blaclaiian 

e) Kaiser 

Out of which Harm and Hainming windows are mostly used 
in many applications. 

3,3 COMPUTATION OF PSD 

There are tiiree methods available to evaluate 

the PSD 

a) To compute the autocorrelation function and then 
finding the Fourier transform. 

b) Direct spectra using Past Fourier Transform 
Technique. 

c) Filtering the signal through a large number of 
band filters and dividing the output through 
each band by the Bandwidth. 

Out of these the second method which finds the 
PSD using Fast Fourier Transform Technique is found to 
be the fastest. Hence the same will be used in this 
application. For smoothing out tiie raw power spectral 
density, Hann window function will be used as it is 
easier to implement on a micro computer. The Hann vjindow 
function in the frequency domain is given by 

Uj (f) = ^ U^ (f ) + ^ Cf - ^ 

C3'*23) 

where U^ (f ) is the modified function 
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(f) is the raw function 
2T is the record length (from, -T iso T) 

The computation procedure for PSD is then will be as 
follows s 

a) Obtain a data sequence of N samples* 

b) Detrend the data using appropriate filter 
(removal of Mean) 

c) Compute the Discrete Fourier Transform of the 
sequence and find out the real and imaginary 
Fourier coefficients. 

N-1 -3211 hi/N 

f Ck) = I X. e = A (k) + a B (k) 

i=0 ^ 

d) Compute the one sided raw power spectral den- 
sity using the following equation 

S,. (k) = I ["A^Ck) •+ 

e ) Using appropriate spectral window on raw esti- 
mates, compute the smoothened Power Spectral 
Density. 

f) Calculate snioothened PSD for each realization 
of the ensemble and take an average to get the 
PSD of track irregularity. 

3.4 FAST FOURIER TRAHSFORM CFFT) Algorithm ril, [2] 

The power spectral measurement involves the 
computation of Discrete Fourier Transform of the data 
sequence. The direct evaluation of the DFT involves N 
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complex multiplications and (N - 1 ) complex additions 
for each value of X (k) and since there are N such values 
to be determined, multiplications and N (N - 1 ) addi- 
tions are necessary. For larger W, direct evaluation of 
DFT involves considerable amount of computation • Fast 
Fourier Treaisform (FFT) algorithm evaluates the DFT of a 
sequence in a efficient way by reducing the total number 
of computations. The FFT algorithm does this by splitting 
the N point sequence into two N/2 point sequences* 

Let the N point sequence be x (n) n = 0, l,,,,!!- 
where N is a power of two. Defining the two N/2 point 
sequences x^ (n) , (n) as the even and odd members of 

X (n) respectively 

x^ (n) = X (2n) n = 0. i , , ~ - 1 

^ ’ ’ 2 (3.24) 

Xp Cn) = X (2n + 1) n = 0,1, ... j - 1 

Then the N - point DFT of .^x (n)^ can be expressed as [I Js' 

X (k) = (k) + ^ X^ (k) (3*25) 

where X.^ (k) and X^ (k) are the N/2 point DFT* s of 
X.J (n) and x^ Cn) respectively. Since X (k) is defined 
for 0 4 k^ K - 1 and X.j (k) , X^ (k) are defined for 
0 ^ k< ^ - 1» To interpret for the values of k , 
the following equation can be used 
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This can be obtained using the periodicity property of 
DFT i.e. = ~W^ • . In a similar manner the N/2 

point DFT can be expressed as a combination of N/ij- point 
DFT, i.e., 

(h) = A (h) + ® 

Ck) = A (k) + ¥|^ B- (k) <^3*27) 

where A (k) and B (k) are the N/H- point DFT»s of the 
even and odd members of x^ (n) respectively* Ttie above 
process can be continued until a tv;o point DFT arrives* 

A two point DFT, F (k), k = 0, i ma/ be evaluated using 
no multiplications as 


F (0) = f (0) + f (1) W 
F (1) = f (0) + f (1) W 


( 3 * 28 ) 


where f (n), n = 0, 1 is the two point transform being 
transformed. 


Since Wg = 1 


if 

and Wg = -1 there are no 


multiplications involved in computing the above equation. 
The total number of operations using this scheme is of 
the order of N log^ N. Thus it can be seen that when 
N is large, there is a great reduction in terms of compu- 
tation. 


The above described algorithm is also called 
as Decimation-in-Time (DIT) algorithm. Since at each 
stage of the process the input sequence is splittedinto 
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smaller sequences for processing, Unother popular form of 
the algorithm is Decimation-in-Frequency (DIF) algo- 
rithm, For this version, the input sequence x (n) is 
partitioned i nto two sequences each of length N/ 2 » The 
first sequence (n) consists of the first two H/2 points 
of X (n) where as the second sequence x^ (n) consists 
of the last N /2 points of x (n) • The same mathematical 
analysis can also he done for DIF algorithm The dia- 

graan in Fig, 3,1 shows the DIT and DIF algorithms for N = 8, 

From these illustrations, it can he observed that, 
for DIT scheme, for the output to he in natural order the 
input should he shuffled where as reverse is the case for 
DIF scheme. Hence whether i.t is a DIT or DIF shuffling of 
the data is necessary. Although shuffled, the order of the 
shuffled sequence can he determined in a simple manner* When 
N is a power of 2, the shuffled sequence will he the hit 
reversed number of the index \ihen it* is represented in Binary. 
For example the following Table shows the shuffled sequence 
for H = 8, 


Table 3*1 


Index 

Binary 

Representation 

Bit Reversed 
Binary 

Bit Reversed 
Index 

0 

000 

000 

0 

1 

001 

100 

if 

2 

010 

010 

2 

3 

011 

110 

6 

if 

too 

001 

1 


101 

: ■ lOt f 

5 

6 

iio 

011 

3 

7 

111 

. ill 

7 
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The same procedure can be adopted to any index 
when N is a power of 2* To calculate the FFT coefficients 
, either the coefficients can be formed as a table and 
then be referenced at various stages during the computation 
or these coefficients can be computed using 

=: Cos ^ j Sin |^( 

where in the first method the extra storage is needed while 
in the second case computation time is involved. 

The FFT algorithms can be further categorized into 

a) Inplace Algorithm 

b) Natural Input-Output Algorithm 

a) Inplace Algorithm; In place Algorithm is the one in ' 
which a given component of any intermediate vector may 
be stored in the same location occupied by the corres- 
ponding component of the preceding vector. This requires 
less storage but computational time is higher than the 
Natural Input Output algorithm. Other characters tic is 
that either the output spectrum appears in an unnatural 
order or that the input data should be arranged before 
entering the computational array, as mentioned earlier. 
This shuffling process is also called as scrambling 
operation, 

b) Natural Input-Output Algorithm; This is the one in which 
a given component of any intermediate vector may not be 
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stored in the same location occupied by the correspond- 
ing component of the preceding vector, thus requiring 
the extra memorjr for storing the Intermediate Results. 
An h-point FFT (H real and N imaginary) will require 
words of memor3^ as compared to 2 N vrords required for 
Inplace algorithm. These algorithms maintain a natural 
input-output order and do not require scrambling opera- 
tion and as such are faster when compared to the inplace 
algorithm. 

Considering the present application, inplace algo- 
rithm is adopted to reduce the memory constraints of the 
system. 

Basic Properties of Inplace Algorithm: 

From the signal flowgraph . the following points 
can be observed* 

a) There are r = log 2 ^ number of computational 
arrays. 

b) In each array every node has two incoming paths and 
two outgoing paths. 

c) In each array,- there are two input nodes whose 
input paths originate from the same pair of nodes 
in tile previous array. Two such nodes are grouped 
as »a dual node pair*. In any array there are 
N/2 pairs of dual nodes* 



56 


d) Each array requires N/2 complex multiplications and 
N complex additions * Therefore the total member of 
multiplications and additions required are 

iTi/z) logg N and N log^ N respectively. 

e) Computation of dual node pair requires only one 

multiplication and tv;o additions. If the weight- 
ing factor at one of the nodes in a dual pair is 
W^j then the weighting factor at the other node of 
the pair is Then 

(k) = (k) + X^__^ (k + N/2^}. 

(k + N/2^) = (k) - (k + N/2^) (3.29) 

"til 

where (k) indicates k component in the 1 
array. 

f ) The spacing between dual node pair diffeis from array to 

th 

array. In the 1 array (1 = 1, 2, • r) the 

spacing is N/2^ means X (k) and X (k + N/ 2^) cons- 
titute a dual node pair* 

g) To evaluate the value of p which is the exponent 

of W for any Index in a given array, represent 

isli 

k, the node index in the 1 array in binary form 
with r bits, retain the most significant 1 
bits and add (r - 1) leading zeros to form a r 
bit binary number. Reverse the bit order of the 
resulting number and the decimal equivalent of 
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the final binary number gives the index p* Ihe 
weighting factor for the k node of the 1 array 
is 

h) The output after r arrays is in scrambled form. 

To unscramble the output X (k) represent the index 
k in binary form with r bits and reverse the bit 
order. Tlie resulting decimal number is the index 
n of X (n), 

3. If IMPLEMEiWATION OF FFT ON 6800 MICROPROCESSOR 

As it was found |^33.» considering the accuracy 
point of view, a 16 bit word length should be chosen to 
represent each data word. For a real time application, 
though less accurate, fixed point arithmetic will be faster 
than floating point arithmetic. Hence fixed point arithme- 
tic is chosen for this application where the data is repre- 
sented in Two's complement form, and the binary point is 
assumed to lie to the- left of leftmost magnitude bit. As 
the execution of the program moves from stage to stage, the 
magnitudes of the numbers in the sequence generally increase 
hence there is a possibility of overflow during different 
stages of computation. To prevent this appropriate tech- 
niques of scaling should be adopted in fixed point arithme- 
tic. 

Scaling |^3^ i The power of 2 algorithm, as mentioned ear- 
lier, operates on two complex numbers. It takes 
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these two numbers and produces two new complesi numbers which 
replace the original ones in the sequence. 

Let (i) and ^3)- he the original complex number. 

The new complex pair is given b^?- 


Vi 


^m+1 


\ hi (3^ 




C3.3O) 


The algorithm goes through the entire sequence of N num- 

M 

bersj taking two at a time. Vftien N = 2 , total number of 

computational stages is M. With the assumption of binary 


point lying towards extreme left, the relation between the 
numbers in m and m + 1 stage is as shown in Fig. 3*2* The 
outside square gives the possible values of Re (i)3 < 1 

and Im [x^ ^iQ 'I the circle inscribed in the square 

gives the region X^^ (i) 1. Similarly the inside circle 

gives the region J X^ (i)| < I/2. If the numbers X^^ (i), 

X^ (j) lie in the region | X^ (i)i < l/2j then X^_^^ (i), 

(0) will be in the larger circle hence overflow condi~ 
tion doesn't arise. However, if X^^^ (i) and X^ (j ) are inside 
the smaller square, it is possible that X^_^^ (i), O') 

to be outside the larger square hence overflo^j arises. So 
the input sequence should be controlled such thatjx^ (l)|< 1/2 
to prevent overflow • There are three techniques of scaling, 


which can be used for this, namely 

a) Shifting right one bit at every iteration 
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Td) Controlling tlie sequence so that | (i)| ^ l/2 

c) Testing for an overflow* 


Out of these, the first method is the simplest 
and easy to adopt for microprocessor applications and hence 
it will iDe used in this implementation. 


A flow chart, and the corresponding loops in the 
FFT in place program is given in Fig, 3 « 3 » stated earlier 
in the FFT routine, the k component in the 1 array is 
given by 


(k) 


X^^^ (k) + (k +. N/2^) 


(3.31) 

X^ (k + N/2^) = X^^^ (k) - (k + N/2^) 

where is the multiplier and is a complex number, given 

by -0 2TT 

= cos (^-) p - D Sin C ) P 


W 


P 


=r e 


N 


As it can be seen, for generating the weights of ¥, 
sine and cosine values are required^ and since calculation 
of sine and cosine functions is time consuming, these values 
can be generated and stored in the system memory before enter- 
ing into the FFT program, for reducing the execution time. 

For this N/2 values of SINE and N/2 values of COSINE aire 
required from 0 to x ( ^ - l) in a stepsize of , 

which corresponds to. N/2 weights of W from VJ®, to 

N/?~1 

W , This requires extra N words for storing the values 


of SINE and COSINE. 
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However, values of sine and cosine can 136 derived 
dy storing only N/4-+1 values of sine from 0 toTT/2 in a 
stepsize of 2TT/H i*e. from W°, W^, ... to as descri- 

ded delow. For N = 1024, this needs the values of sine to 
de stored from 0 toTT/p in a stepsize of 271/1024 \^hich 
corresponds to v7°, to 

Calculation of Sine Values (N = 1024) 

a) First whether the index ^ 256 or > 256 is to de checked, 
d) If index is ^ 256, then use the index as it is to fetch 
the appropriate sine value from the array. 

c) If index is >256, -the index is subtracted from 5l4 (in 
general from 2 (N/4+l)) and this value is used as index 
to fetch the a,pproprtate sine value. 

Calculation of Cosine Value (N = 1024) 

a) The program should check whether index ^ 256 or y 256 
(in general index ^ or ^ ). 

d) If index is ^ 256, subtract the index from 256 (generally 
from N/4) and this value will de used as index to fetch 
the value from Sine array. 

c) If the index is > 256, subtract 256 (generally N/4) ^ from 
the index and this value will be used as tlie index to 
fetch the value from Sine array which gives appropriate 
values. 
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By this methodj the memory required will be 
N/if+1 words for storing Sine values to generate the weights 
which will cut down the memory from N words to NA-+1 words* 

Organization of the FFT Program 

Bae FFT subroutine assembly language program mainly 
calls the following subroutines; 

a) COMPmiEll!: ROUTINE (COMPL) 

This routine finds the, Two *s complement of 16 bit 
number when called* 

b) OVERFLOW ROUTIlffi (OFL) 

This routine shifts the 16 bit number to right by 
one bit, equivalently it divides the number by 2* 

c) POWER OF TWO ROUTINE (LOGpN) 

This finds the log„ N where H is the number of 
samples and is a power of Two* 

d) DOUBLE PRECISION MULTIPLY ROUTINE (DPMUL) 

This multiplies the two 16 bit numbers and produces 
a 16 bit product using Booth’s algorithm. Negative 
numbers are represented in Two’s complement form, 

e) BIT REVERSAL ROUTINE (BITR) 

This finds the Bit reversed number of a 16 bit 
number • 

f) UNSCRAMBLING ROUTINE (UNSCM) 

For inplace algorithm, since the output is in 
scrambled form, this routine orders the output 


sequence • 
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A complete assembly listing of FFT for 6800 micro- 
processor applications v;hich can be used upto N = 102^ is 
given at appendix. The listing to find the power spectral 
density which calls the FFT routine is also included in the 
Appendix. The memory requirements and the results are pre- 
sentedj details are discussed in Chapter 5* ■ 
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MiCROCOi'ffUIER BASED SYSTEM 
4.1 INTRODUCTION 

The technological revolution, spearheaded, by the 
MICBORROCESSOR, a single chip central processing unit has 
rendered cost effective real time signal processing practical. 
A typical real time signal processing scheme like the track 
profiles measurement using large main frame computers or 
minicomputers would have been an expensive proposition and 
therefore impractical for routine track evaluation. With 
today’s microprocessors, a low cost, completely independent 
on line dedicated signal processing system can be built for 
such applications. 

The 6800 microprocessor and its family support is 
chosen for the instrumentation proposed for the real time 
evaluation of railway track profiles. In this Chapter, the 
track profile measurement system’s specifications, system 
architecture, and hardware details of the system which also 
include, system packaging have been described in detail. 


4,2 SYSTEM SPECIFICATIONS 


To make the system upward compatible and to have a 
proper system-to-user interface and also for having electrical 
compatability with the existing rail cars , the following 
specifications fdr-' tfee--y^|.'';’%sfeK3a6niihrifig‘"''sysiem have been 
laid down. ■' '■ fV::’ , , - ■' 
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4.2.1 Environmental Specifications 

a. The system should he able to operate on + 24V DG battery 
system normally available on railway carriages. 

b. It should be portable and should not require any air- 
conditioning . 

c. It should be able to survive in the vibration environment 
normally obtainable on a coach floor. 

4.2.2 Signal Processing Requirements and 
Real Time Nature of the Signal 
Processor 


a) Track Parameters to be extracted: For every kilometer 
span of track length, absolute vertical profiles of 
left and right rails, absolute lateral profiles of left 
and right rails, dynamic cross level, centre line uneven- 
ness, dynamic gauge, centre line alighment and power 
spectral measurement of each of the track profile is to' 
be evaluated by the system in real time (i.e. max. 3^ 
seconds : the time talcen for the rail car to traverse 
1 Km at the envisaged maximum speed of 120 Km per hour). 


b) Analog Signal Channels: The system is to process the 

outputs of four accelerometers mounted on the axle boxes 



for both horizontal and lateral profiles of the left and 
right rails. Hence tho system should be able to sample, 
digitize and ..from -these four.i;; 
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c) Conditioning of the Input Signals? The outputs of the 
accelerometer should he suitably conditioned to the 
compatible levels of the system input digitizer. 

d) Number of Points Per Channel? The number of sample points 
for each channel is fixed to be 512 for every kilometre 
section of the track vmich corresponds to an approximate 
spatial samp3.ing of 2 metres. These sampling instant 
pulses are generated by the wheel and its associated 
electronics. 

e) Processing Time: Considering the maximuii speed of rail 
car to be 120- Kilometres per hour, there is 30 seconds 
of time available for each Kilometre traverse. Hence 
for each Kilometre of track span, all the four channel 
signal processing and track categorization should be 
carried out within 3 0 seconds. 

f) Machine to Operator Interface? The system should be 
capable of accepting the Initial kilometre entry at the 
time of starting of the run 9,nd automatic updating of the 
kilometre should be done on Command at the start of each 
kilometre of track. A.utomatic recording or printing of 
the results at the end of tKe run, as may be required, 

is also nece 

g ) Track to 
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conditioned before feeding to the signal processor. 

If. 3 SYSTEM DESCRIPTION 

Tile system architecture is shovm in Fig. i-f.i. It 
comprises of a general microcomputer system and a monitor 
on a bus and application dependent modules. It consists of 
the following module si 


a) 6800 processor based CPU module with monitor ROM, stack 
memory, serial I/O (ACIA), parallel I/O (PIA). 

b) l6 K ROM module for program storage. 

c) 16 K Dynamic RAM module. 

d) Ram Refresher Module which generates all the timing 
pulses for Dynamic RAM module and address decoding 
circuitry. 

e) 8K Static RAM module (optional ) 

f ) Digital Hardware multiplier module. 

g ) Console Module which consists of a fifteen digit eight 
segment display and a Hex Keyboard. 

h) Cassette and RS 232 C interface module. 

i) 20mA Loop TTY interface and Baud rate clock module. 

5 ) Floppy disc drive controller module ■ (optional ) . 

k) DM controller module (optional). 

l) Analog output device’ controller module (optional). 

m) A/D converter mpauiB-. '■'iY; ; ^ ^ 


n ) Analog =7; 








4 ^ 
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In the above, the modules Ca) to (k) are for 
configuring a general microcomputer system ^toore as the 
modules (l) to (n) are for the envisaged specific applica- 
tion. 

Microcomputer System Hardware 

In this section, the modules which comprise the 
microconiputer system that will be used in the proposed system 
have been described. To make the system modular, double 
sided printed circuited boards of size 6 - 3 -" x have been 
used for all the modules v;ith a 44 pin edge connector con- 
tacts that go on to the system bus, Ihe edge connector 
tonguG designation with associated bus signals are given in 
Table 4.1 and the same will be followed to all the hardware 
mDu.«.Xoc! except whei'e mentioned* 

4. 3 . 1.1 CPU Module 

The complete circuit diagram of the CPU module is 
shown in Pig. 4,3* This mainly consists of a CPU clock and 
decode circuitry, on board baud rate clock, processor, memory 
and I/O provisions, 

a) CPU Clock and Decode Circuitrys Since the system uses 
M6800 microprocessor which uses a two phase clock 
(^l, 02 ^ i 44 is required to generate a two phase non 
overlapping clock that should run at V level, accord- 
ing to the requirement of, 6800 processor., ,The M6800 
uses this' a./4 MHz clock has ' 7 :7'.. ',v/ 
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iDeen used comprising of a h Mz crystal and a part of 
U1 (7^00 IC), tile output of which, v/ill bo further divided 
by a factor of four using U2 ('7173). 7he output of U2 
will be a 1 Mz clock which will be buffered using 
U3 (7104-), part of U1, U5 (7437) and two transistors 
^2 ^2W2222) to derive a nonoverlapping two phase 
clock. Tne decoding circuitry consists of U4 (74 i 0) 
wliich generates tiie various timing signals, .using Valid 
Memory Address (VMA) of the CPIJ, like VNA'.J^^ j 

RCS , The RCS is the chip select for the Read Only 
Memory (ROM) on CPU module where the 'Monitor' of the 
system actually resides. RCS is 'AMDED' configuration 
of 02 clock, R0MEN1, R0I-1EN2 where R0MEM1 is ROM ENABLE1 
and R0MEH2 is ROM ENABLE2* These signals are derived 
from Ram Refresher Module. ROMENi is a buffered version 
of R/Vf (Read/write) and R0MEN2 will be in active state 
when the MSB address lines have got an Hex address of 
'F'. 

b ) On Board Baud Rate Clock.: This clock will be used in 
conjunction i/itii Asynchronous Communication Interface 
Adapter (ACIA 6850) for serial I/O application. This 
is a simple astable multivibrator using U15 (555)? a 
timer IC with a precise baud rate- (3OO baud) adjustment 
using a trimpot. If higher baud rates are required 
this on board baud rate clock should be disabled, and 
external clock should be cikinecte.d to the system. 



7 ^ 

c) Processor, Memory and I/O; This section is the heart of 
the s.ystem consisting of 6800 microprocessor, 1 K R/v,f 
memory (a part will he used as stack memory), two Periphe- 
ral Interface Adapters (PIA 6820) for parallel I/o commu- 
nication, an ACIA (6850) for serial I/o apiDlications and 
a 1 K ROM space where the 'Monitor’ ROM resides. The 
processor address, data and control lines are brought on 
to the edge connector contacts which go into the system 
bus. Hie PIA outputs are available on a 50 pin Eurocard 
type connector. The ACIA connections are brought on to 
a "l^- pin DIP socket. For serial and Parallel I/O communi- 
cations proper jumper assemblies should be used for the 
Eurocard type connector and DIP socket. In the system 
configuration that is used for this application, the DIP 
jumper assembly goes to the RS 232 and cassette interface 
module which in turn can be used for the Asynclironous 
communications like cassette, RS 232 C, Teletype etc. 

The : memory mapping of the vai'ious I/O extensions and 
the memory on GPU module is, shown in Table 4*2. The sys- 
tem reset is hard wired through a puslibutton switch to 
the reset pin (pin 40) of the processor (or tongue Z of 
the edge connector).- The serial connector (l4 pin DIP 
socket) signal surrmiary is given at Table 4.3* After the 
power-on condition and the system reset, the monitor that 
resides in the ROM takes .the control of the system. For 
its use, the monitor mainly uses the stack area on board. 
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4. 3*1*2 16 K ROM Module 

This module proTides 16 K byte memory accessory 
for the control module. This laodule uses 2708 erasable/ 
programmable type ROM ICs. The ROM card may be configured 
into any qusulrant of the processor's address space using jum- 
per selection of the higher order bits. The ROM array is 
fully buffered from fne microprocessor bus using low power 
buffer gates to assure minimum loading. The complete cir- 
cuit diagram is given at Fig. 4- .4. ' Ul and U2 are low power 
gates which isolate the capacitive loads of the ROM array 
from the S 3 ''stem bus. to (the LSB address) are supplied 
directly to the 2708 address loads. A^q’ ^12 ^13 

are buffered by ^nd U]^ and used as inputs to a 74-154- (one- 
of-l6 decoder). This decoder generates the chip select sig- 
nals to select one of the 2708 's from the array. Address 
and A^^ are decoded to enable the data bus drivers’ "tJ6 
and.U7* The state of A^]_^ and A^^ which enables the board 
is jumper selectable. Tiie placement of these jumpers will 
determine v/here tho ROM board appears in the processors memory 
map. Four jumpers arc used to malce the memory space selec- 
tion. Only two of 'ciiese juiiipers are ever in place at any 
one time. Jumper insertion and the resulting memory map 


is shorn below: 
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Meiaory locations Jumpers in 

occupied place 


COCO - 


J1, 

J4 

8000 - 

BFPF, ^ 

15 

J2, 

<73 

4000 - 

7FidM . 

1 0 

J2* 

J4 

0000 ~ 

3EPF^6 

J2, 

P3 


4*3. 1.3 I6K Dynamic RAJ'-I Module 

This RAM module is a high density momory accessary 
for the .CPU module. The complete circuit diagram' is shorn 
in Dig. This RAM module uses 4096 dynamic RAM IC, This 

module accepts the multiplexed addresses and timing signals 
from the RAM REIUIESHER Module through the refresh "bus and 
communicates data thru system bus. Each ^096 byte bank of 
RAM has an associated TTL gate to supply the ROW Address 
Strobe (RAS) and Chip Select CcS) signals to the eight IC’s 
comprising that bank. The RAS decoding is used to minimise 
the quiescent power dissipation. This is- achieved by giving 
the RAS signal only to the particular RAI''I bank selected. 
During tlie refresh cycle the REF signal from the refresher 
module forces RAS to be applied to all the memory ICs. During 
the refresh cycle, processor does not receive hence does 
not generate This signal is gated through CS to the 

RAMs so that all RAMs are deselected during a refresh cycle. 
When RAI4 refresher module is used, the onboard clock shoiild 
be disabled by removing U2 on GPU module since the clock is 
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now supplied by refresher raodulo. Half of is used to 
detect a reference to the memory board. Pin 6 of this IC 
goes high when any of the four banlvs is to be read from or 
i^ritten into. The other half of U4 and part of U7 generate 
enable signals for a pair of Tri-state transceivers U5 and 
U 6 , which have iiigii input impedance to minimiae the loading 
on memories and processor data bus. Only one half of these 
transceivers is enabled at eaiy one time. The half that is 
to be enabled is determined by the state of the R/W line from 
the refresher module. This R/W signal is also combined with 
the WRITE pulse from the refresher module to apply the ^^rite 
signal to the RAM array at the proper time in the memory 
timing cycle. Most of the U 3 is used as a buffer to drive 
the capacitive load of the Column Address Stroho (CAS) 
pins on the RAM array whore as U2 is used as buffer for the 
multiplexed addresses for the RAM array. Typical wave forms 
of the timj.ng signals are shoi/m in Fig. 4.6. These signals 
are buffered by TTL gates to minimize loading on the signals. 
Tile configuration of the RAM module in the processor ’s 
memory space is jumper selectable. Each 4096 bite 
bank may be configured independently at sixteen locations. 

The jumper termination with its associated address is given 
in Fig, 4 . 7 . All the supply voitages are derived from 
the Ram refresher module, . 
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Rara Refresher Module 

This module is the linl-i "between the Dynamic RAM 
module and the CPU module. This provides the required timing 
and refreshing signeRs for dynamic RAM module and also "buffers 
the address signals from CPU module. The complete circuit 
diagram is shown in Fig, ^.S. 

IC1 comprises an 8 MTz crystal oscillator for the 
processor master clock which will "be o"btained "by dividing the 
8 Mz clock (pin 8 of IC1). When this module is to "be used 
with CPU modu'ie the clock on CPU module should "be disabled by 
removing the UIv flip flop (IC2) on that module, 

IC5 and IC6 (74-95) form a divide-by- eight counter 
and meraory clock generator. The signals shown in Fig, 4-. 9 
with their X'jraveforms are generated by this shift Register 
chain. The Row Address Strobe ( RAS ), Coulmn Address Strobe 

( CAS ) and WRTTE are the required timing signals for the 
dynamic RAM module (refer description of Dynamic RAM module), 
the processor clock is supplied to the system bus through 
a Nand Gate (IC7) and is also used to operate the memory 
address multiplexer consisting of IC10 and IC13« These ICs 
are tri-state low power schottky devices used for minimizing 
the bus loading. IC2, IC3 and half of IC4 form a 32p.s 
timer which determines the refresh interval. When a refresh 
cycle is required, pin 12 of IC4- goes high; 125 ns into the 
next memory cycle pin 9 of IC4- goes high indicating that 
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the memory^ cycle is to be a refresh cycle. Refresh addresses 
from IC85 ^^9 are supplied to the memory board through IC11 
and IC12. The clock is suppressed during the refresh 
cycle. The rlER signal is given to the memory card to signal 
it that a refresh cycle is in progress. The higher order 
bits of the address -^^3? ^^1^5 ai'’© buffered by 

IC17 and are presented to a lour-to-sixteen line decoder 
(ICl4 7^15^). The output of this decoder is supplied to the 
RAM card where these will be used to select 4 K byte RAi4 
banks. IC14- asserts 'RA14SEL when it decodes 'E' on the 4 
highest address bits. The CPU module and interface modules 
use RAMSEL along with decoded address bits to generate chip 
select signals to activate the various PIAsj ACIAs etc. U15 
asserts R0iiEN2 when it decodes ‘F» ©n the four highest address 
bits. R0MEN1 is the buffered version of R/W. Both R0MEN1 
and ROMEN2 are used to enable the ROM on CPU module. 

8 K Static Ram Module 

Although packing density is low, the advantage with 
static RAMs is that refreshing is not necessary for these 
memories. - A 8 K static RAI4 module using Intel 2114 (l K x 
4 bit) static RAM is also designed and the complete circuit 
diagram is shov/n in Fig. 4.10, In . this all the address 
lines are buffered using low power schottky devices (to 
minimize bus loading) IC3} IC4, IC5 (74LS365) •R/VJ, ' V14A."^ 
are also buffered using these ICs. Data lines go to two 



80 


data bus traiiscGi'vers (8833) tliat coi-rmiunicate v/itli meraory ICs^ 
In tliiSj there are two 4 K memory banlis and the memory space 
allocation for these banks is jumper selectablic, Hiis is 
done using two 3'“to~8 line decoders (7^155) and two OR gates » 
and the outputs of 'diese OR gates will form as chip select 
signal for a i K x 8 bit section in which two 211^s will be 
paralleled. The enabling of the data bus transceiver is done 
using the two bank select signals and R/V/ signal; R/W selects 
the mode of data bus transceiver i.o, transmit or receive. 

This module is also bus compatible following same tongue 
designation for the edge connector contacts. 

!+,3,1,6 Console Module 

llie console I/O module provides an economic means of 
directly inputting / outputting information between operator 
and system bus. The complete circuit diagram is given at 
Fig. 4-. 11. The module has got two main parts, display and 
keyboard. It has got a real time clock which interrupts the 
processor at regular intervals. The display consists of 
fifteen digits of eight segment RED readouts and the display 
codes are mider full sof tv-rare control for maximum versatality. 
The refreshing feature of the display eliminates the need 
for latches and accomplishes this by turning on one digit at 
a time. The segments that are to be turned on will be de- 
termined by Ill’s side” bits. These are buffered 
by U5 and U6 wwich then drive the segments through current 
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limiting resisters through Rg. Each common cathode of 
the readouts is driven hy U3 and U^i- , one-of-l6 coder 
driven 'hj the display refresher counter U7 (7^93)* U7 is 
cleared by CA^ of U1 . No digit is on while U 7 is clear. 

Then each successive digit is turned on, from left to right, 
by strobing U 7 »s clock input by CB^ of Ul . Control register 
B of U1 is therefore configured so that when a new set of 
data is written into the PIA, GB^ goos low in order to switch 
on the next column. The display can be turned off completely 
by control of CA^ of U2. A ’zero' turns on the segment. 

The PIA (Ul) A-side with its internal pull up 
resistors is connected directly to the key switch. The 
common row terminals of each switch are connected to PAh-PA 7 
of Ul, The common-column terminals are connected to PAO-PA3 
of U1. If the A side is configured as an input, the key 
switch may be inputted by reading the eight bit A-side data. 
U8 (CD^ 060 ) is a real time clock (CMOS oscillator /ripple 
coimter) which gives interrupts to the CPU at a rate of 75 > 
150, 300 or 1200 interrupts/second depending upon jumper J3 
connection. Normally IRQ Interrupts are generated by the 
console module. However a simple jumper modification allows 
ilMI type interrupt. The display refresh routine will be as 
follows. The routine receives interrupts from the module 
configured for an interrupt interval consists of two sec- 
tions; Initialisation code performed once at reset time 



and inter3?up-u drj.ven code which at interrupt time 
display code foi* the next digit and outputs it to 
drirers . Sie memory space occupied by the module 
■via jumpers 11 and J2. 


gets the 
the display 
is S’clocted 


4. 3* 1*7 Cassette and RS 232 C Interface Module 

This module pro-vidos tlip’ asynchronous communications 
to the GPU module through the use of a cassette deck or RS 
232 console. The complete circuit diagrara of ‘bhis raodule is 
sho'v7n in Fig. 4-. 12. This module is 'no'i bus compatible' 
hence should bo mounted separately, and uses the.AGIA jumper 
assernbl;,^ frem CPU module. 

WRITS Operations During ca.ssotto write operations, the RTS 
lead froia ACLi is made low. This makes the IC14 as a dividc- 
by~four or a divido-by-eight counter. The divide constants 
depends upon 'bh^j s’Cate of th.; TXD lead from AC lA. The 4800 Hz 
clock (300 baud clock) from casset'be clock and TTY interface 
board is divicud to 2400 Hz for a 'one ' from the ACIA or 
1200 Hz for a 'zero B'or high speed operation the 2400 Hz 
ou'bput of IC5 is supplied to the ACIA operating in a divided 
by 'one' mode. Lov; speed operation is achieved by s. electing 
the 4800 Hz clock while operating the ACIA hi -011.0 divided 
by sixteen mode. The timing chart for m.‘±tCi operation is 
shown in Fig, 4,13* 

As it can be seen froia the timing chart for the 
high speed operation, 'OIE ’ is -recorded as a single cycle 
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of 2^C’0 Hz v/,aile ’zero' is recorded as one half cycle of 
■1200 Hz» In tliG low speed operation ’one’ is recorded as 
eight cycles of 24-00 Hz while zero is recorded as four C3^cles 
of 1200 Hz. Cassette write operations inhibit the output to 
the RXD lead by allowing IC5 to be clocked set thus forcing 
RXD to’ the ACIA to be held marking* Output to the RS232 
console is also inhibited by IC 8 * 

READ Operation! IC2 and most of IC 1 and IC 3 comprise a 
limiter /edge detector v/hich squares the tape audio data and 
generates short positi-ve pulses at the output of IC 3 (pin 11 ) 
which is also the reset signal to the baud rate clock. Hence 
these pulses are used to synchronize baud rate clock to the 
incoming data. The clock signal (300 baud) from baud rate 
clock is supplied to a missing pulse detector comprised of 
IC4. Missing pulses indicating ’zero’ data force the output 
of IC9 (pin 11 ) high which causes the R]Q} lead low for each 
zero bit detected. Read clock is recovered from the data by 
IC5. Phasing ambiguity is removed when reception of the first 
zero occurs. The missing clock pulse allows a low to be 
loaded into IC4 (pin 9) • IC9 pin ^ is constrained to become 
high in the next transitioh of the recovered clock so that 

■ , i ' 

the succeeding bits will be centre sampled by ACIA. Input 
from thf: RS232 console is locked out by the data multiplexer 
comprised of IC 6 and IC9* Output to the console is inhi- 
bited by ICS. Read timing signals are shown in Fig. 4.l4. 
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Tnis cassette and HS232 interface module is 
also provided with a logic to changeover the connections 
between RS232 and TTX. Since TTY uses tiio same ASIA jumper 
assembly this logic is necessary which comprises of IC9) 
ICIO5 IC11, IC125 ICI3 as shorn in the diagram. The RS232 
interface comprised of IC8, IC15 (MCl!+88 Imie driver) and 
a transistor (2N 2222) and D^. 

^«3»1*8 Baud Rate Clock and TTY Interface (20mA Loop) 

This module provides the necessary baud rate clock 
for asynchronous communications and also has got a TTY 
interface (20mA loop) circuitry. Hie circuit diagram is 
shoi'm in Fig. 4,1 5 and this module also not bus compatible. 
Hence should be mounted separately. 

Baud Rate Clock; This comprises of a crystal controlled 
307.2 Kliz clock derived 13'" two inverters (IC5) and three 
delay type flip flops (ICl, IC2) vmich will be used by 
cassette interface module. This master clock is further 
divided by 103? IC4 vriiich form a. counter and give standard 
terminal baud ra/bes viz. 150} 3OO5 600, 1200, 2400, 4800 and 
9600, These baud rates from 150 to 9600 baud are switch 
selectable using a DIP type switch. These switch selected 
rates are supplied to the AGIA via a clock multiplexer cir- 
cuit comprised of IC6 and IG9 on cassette interface module. 
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TTI Interface (pO mA loop) s The necessarj^ TXD, RXD signals 
for serial transmission of this TTY interface are derived 
from a. simple switch over logic between RS 232 -' TTY on cassette 
interface module. The baud rate will be fixed at 3 OO baud 
by connecting the on board baud rate clock on CPU module to 
the TXC, RXC clocks. Tliis is also done by the switch over 
logic* The T](D signal is buffered by a UAIID gate (IC 6 ) and 
is coupled to a opto isolator (IC7) which in turn supplies the 
necessary current for the TTY loop. Similarly TTY output 
signals are coupled through a opto isolator (ICS) and a MMD 
gate CIC 6 ) to the RXD terminal of ACIA, The TTC. logic levels 
are converted to a current loop by the opto isolators. 

if, 3 , 2 I'licrocomputer System Configuration and Software 


The above described modules configure a complete 
microcomputer system with the input/ output provisions for a 
cassette deck, RS 232 console or a Teletype and a console 
module with a Hex key board and fifteen eight segment LED 
display. The basic system memory map will be as follows; 


a) CPU Module 

Hex address 
0000 - 03PP 
ED80 - EDPF 
EE08 
EE 09 


Component 

RAM 

RAM: (stack .memory) 

U9 ACIA control/ status r eg is tei* 
U9 Data Register 
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EE10 

U8 

(PIA-1 ) output register A 

EE1 1 

,U8 

output register B 

EE12 

U8 

control register A 

EE13 

U8 

control register B 

EE20 

U7 

(PIA- 2 ) output register . 

EE21 

U7 

output register B 

EE22 

U7 

control register A 

EE23 

U7 

control register B 

FCOO - FF] 

FJ' Read only memory 


Td) Console module base address at EESO^^^ 

c) 16K I^mamic Rain memory allocation from #MNI 

d) 8K Static Ram memory allocation frora I®#* * If f P '111 

e) 16K ROM module memoiy allocation from — 

f) A/D converter PIA address 

Software 

flic system software supports a ROM based lK-byte 
loader/monitor/debugger programme for 6800, Tlie monitor 
resides to 2708 ROM and is on CFJ module. The monitor supports 
the following coi;amands toen entered from a teletj^^pe or video 
terminal. 

Features of tiie- monitor*. Thee following are the monitor 

commands* 

M - Display/change memory 
I - Display/ change instruction 
- Display processor registers 


R 
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CNTL-A - 
CNTL-B - 
CNTL-C *- 
CHTL-P - 


Displa^^/ change accumulator A 
Display/ change accumulator B 
Display/ change condition code register 
Display/ change program counter 


GNTL-X 

G 

S 

D 

P 

L 


Display/ change index register 
Go to user program 
Single step user program 
Dump memory 

Punch loader compatible tape 
Load memory 


The monitor also supports 

. free-format hexadecimal input 

• insertion and deletion of brealc-points 

. interrupt vectors 


iTee-format means the leading zeros can be omitted 
when entering hexadecimal address, data, OP codes, or the 
first or second byte of an operand. After the power-on 
condition and pushing the reset switch reverts system control 
to monitor. Then the monitor responds with a carriage 
return (CR), line feed (LF) and an asterisk (*). The 
indicates that monitor is in the command mode and ready to 
accept any of the monitor commands. RESET should be used 
whenever power is first applied or system control to be 
returned to monitor. 


Interupt Vectors 


a) Software Intoriipt; If the processor encounters an SWI 
instructioiij monitor responds by pushing the MPU regis- 
ters onto the users stack and then vectors to memory 
address EDFA^^* Memory addresses EDFA, EDEB and EDFC 

are reserved for the users SV/I handler or a jump ins- 
truction to the SWI handler sit another memory location* 

b) Interrupt Request; If the IRQ line is pulled low and 
the interrupt mask bit is not set, monitor responds by 
pushing the MPU registers onto the user stack and then 
vectors to memory address EDF7^^. Memory locations 
EDF7} EDF8j and EDF9 are reserved for ‘fclie user IRQ 
handler, or a jump instruction to the IRQ handler at 
another raemorj;" location* 

if. 3* 3 The Specific Application Dependent Hardware 

For the track profile measurement, since the signals 
ape in ^alog fo 3 :*rii A/D converter is necessary for digiting 
the input signals. Hence the following modules ore used”. 

a) A/D convertor Diodule h 

b) Analog multiplexer module 

Converter Module 

The complete circuit diagram of this module is 
shown in Fig. 4.1 6 whicli is configured using 8700 series 
A/D converter i In this, the eight bits of converted analog 
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data are presented to PIA. peripheral bits PA7 to PAO. The 
most significant bit is PA7 and least significant bit is PAO. 
The A“Side control lines give some freedom of control over 
the ADC. The ADC may ’free run’ or convert only when required 
under program control. An interrupt can be specified at the 
end of the conversion cycle or programmed I/O techniques can 
be used. To sample a specific channel the following routine 
is adopted. The routine is in two parts i initialization 
routine called only once at the start of the program and a 
sampling routine called when a sample is to be taken* 

The memory space occupied by the module is selected 
via jumpers H and >12, System gain may be adjusted using Pot 
P11 on A/D converter module. The adjustment of this Pot 
allows the ADC to accept and convert voltages from a 1 volt 
to 20 volt interval. Offset Pot PlO allows the ADC to accept 
either unipolar or bipolar signals. If bipolar signals are 
chosen the inherent coding is offset binary. If two’s comple*- 
ment coding is needed, the conversion is done in the software. 
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If .If SYSTEM OPERATION 

The signal processor for its operation, mainly 
uses the feature of interrupt vectors of the monitor. 

The console module uorks as a control panel with some keys 
marked as 'TtlT (Rail Track) monitor", "KM Scan", and ' 
"STOP". Tjfter the power on condition and system reset, 
the monitor takes over the control and waits for an input 
action* The console module, as described in the descrip- 
tion of system hardware, gives regular interrupts. After 
the first interrupt is encountered, the monitor vectors to 
the one of interrupt vector address location. Since the 
console module is configured to M4I tyge interrupt th-e 
monitor vectors to the address location EDFD where EDED, 
EDFE, EDFF are reserved for interrupt routine handler or 
a jump instriiction to the IRQ handler. This routine tests 
the key pressed and if it is RAT monitor key, then the 
routine branches to the track monitoring routine, other- 
wise it waits until the key is depressed as mentioned in 
the operator-system operation interface flow diagram 
(Fig, If # 2 )* After the RAT monitor key is depressed, the 
track monitoring routine displays the word "FEED" on the 
console module lED display. Then 'the initial kilometer 
is entered and the routine performs the initialization of 
multiplexer and A/D modules. Data memory is divided into 
two buffer aroasc -'When the' data is collected, sampled 
and being stored at one buffer, ttie date processing on 
the other buffer is simultaneous}# performed. After 1 Km 



track span data is evaluated, the results are huffered 
into the output data memory and cleared. Then next Km data 
is stored in this data buffer, processing of the data 
already stored in the other buffer is performed. Thus, 
the data buffers are alternatively used for data collec- 
tion and processing. 

The complete system operation is as follows. 

After the initialisation of the multiplexer and A/D 
modules, the routine waits for the keyboard action, and 
checks whether Km Scan Key is depressed or not. This 
key is depressed at the start of the run, where the data 
is to be collected for the processing. If the Km Sean 
key is pressed, the routine increments the Km (present 
Km) and enables the data sample pulse generated by the 
rail wheel and its associated electronics. During the 
first ICm of run, only data will be collected and stored 
in buffer 1 since no processing of the data is involved. 
In the next Km, data is stored in buffer 2 while the data 
in buffer 1 is processed for track evaluation purpose* 

The data processing scheme involves in finding the track 
profile of both left and right rails by filtering and 
double integrating the conditioned output of accelerometer. 
Then power spcictral meas^Jirement of the track profile is 
done and the results are stored in the output data buffers, 
and after this the buffers will be cleared. In the third 
Km run, the data is stored in buffer 1 while the data in 



Touffer 2 is TDeing processed* This alternate cycling 
"between the data buffers is continued until the ’’Stop” 
key is depressed. After the depression of stop key, the 
routine disables the sample data pulse from the rail 
wheel, but continues processing of the data for the last 
Km, Then depending upon the status of the print key i.e. 
if Print Key is depressed or not the results will be prin- 
ted on a printer and the control is taken over by the 
track monitoring routine for the next run. If the suspen- 
sion of the activity of the processor is required at any 
moment, the system reset should be used for this purpose. 
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Jumper Termination and Associated 
Address for RAM Module 


Jumper Termination 

Memory 

Space 

0 
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- 

OFEP 

1 
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- 

1FFF 

. 2 

2000 


3FFF 

3 

. 3000 

- 

3FFF 


If 000 

-tori* 

1+FFF 

5 
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- 

5FFF 

6 
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- 

6FFF 
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- 

7FFF 
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“ 
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A 
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mm ' 

AFFF 

B 

BOOO 

pm 
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*• 
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Table Edge Connector Signal Summary 


Component Side 

V/iring 

Side 

Tong^ue 

Designation 

Signal 

Tongue 

Designation 

Signal 

1 

Gnd 

A 

A1 1 

2 

+57 

B 

A12 

3 

D5 

G 

A10 

4 

D? 

D 

A13 

5 

D6 

E 

A9 

6 

d4 

E 

A14 

7 

D3 

H 

A8 

8 

D2 

i 

ill5 

9 

E1 

K 

A? 

10 

DO 

L 

A6 

11 

-12V 

M 

A5 

12 

+12V 

N 

A4 

13 

MM SEDECT 

R 

A3 

14 

- 

R 

A2 

15 

- 

S 

A1 

16 

ROM EHABLE 2 

^';T ■ 

AO 


ROM EMBLE 1 

5u 5'' , 

B.A 

18 ; 

vm • ■ 


R/W 

19 

BUS 7^2 


EM 

; 20 

stretch 

X 


21 

+5V 

Y 

halt' 

22 

Gl®. 

z 

RESET 


TABLE if, 2 


CPU Module Memory Mapping 


Hex Address 


Component 

0000 - O 3 FF 

HAM 


ED80 - EDPF 

RAM 

(Stack Memory) 


EE 08 

U9 ACIA Control/Status Register 

EE 09 

U9 Data Register 

EE10 

U8 (PIA-1 ) Output Register A 

EE11 

U8 Output Register B 

EE12 

U8 Control Register A 

EE 13 

U 8 Control Register B 

EE 20 

U7 (PIA- 2 ) Output Register A 

EE21 

U7 Output Register B 

EE 22 , 

U7 Control Register A 

EE23 

U7 Control Register B 

FCOO - FFFF 

Read only Memory 



Table ^.3 


ACIA Serial Connector Signal Summary 


PIN Number 


Signal 


1 

2 

3 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 - 


+5V 

On board Baud Clock 
Gnd 
Gnd 
Gnd 
+ 12 ^ 

-12V 

Data Carrier Be teet 
Output Data 
Request to send 
Output clock 
Input clock 
Input data 
Clear to send 



CHAPTER 5 


SOFTWAHE IMPLE>1ENTATI0N, RESULTS AMD 
DISCUSSIONS 

5*1 SOFTLfARE IlvIPIBIENTATION 

In tills section, the implementation of composite 
iiighpass-integrating filter, Fast Fourier Transform and 
PSD measurement of track profile for 6800 microprocessor 
applications is described* A complete assembly program 
listing for these functions is given at Appencix ’A«, 

The memory req.uirement for the individual programs is 
given here. For higher flexibility the data memory loca- 
tions for filter and PSD are considered separately for 
testing* To make up the complete track monitoring rou- 
tine, these locations should be changed properly* 

5*1*1 High Pass filter aid Integration 

This program takes approximately 1K of the 
memory ( 114-6 locations) aid is implemented in 32 tii't 
organization for higher accuracy. , This is assembled from 
Memory address 5660^ ^ to 54-77^ g* The data buffer consists 
of four sections and is assembled from 2000^ to 27FF^£^ 
(Buffer 1 ), 2900^g to SQFF^g (Buffer 2), 3260^^ to 39FF.|6 
(Buffer 3), 4-8EF^g ' (Buffer 4-), each a 2K section 

of the memory. The data stored in the buffer 1 is first 
filtered, integrated and then stored in buffer 3 
again the data stored in buffer 3 (first high pass and 
integration) is again processed to get the double '' 



integrated version of the original data and is stored 
in buffer 1* Then the control goes over to Buffer 2 for 
processing the next section of the data. Buffer 2j 
Buffer h are used for this. After double integration of 
the data in buffer 2> the control again goes back to 
buffer 1 and this repeats. 

5*1.2 Bast Fourier Transform 

The FFT routine takes approximately K of 
memory (^580 locations) including the data locatipns for 
real and imaginary parts each a 1K section (for 5l2 points) 
and is assembled from 4-000^^ to 5lFV^gi» The data loca- 
tions for real part start from to 4-DG9^5 and 

4l)CA^^ to 51C9^^ are for imaginary part. The data on 
which FFT analysis is to be done is first stored in these 
locations and after the run, the result is available at 
the same locations* 

5*1 #3 PSD Measurement 

This routine takes 353 locations of the memory 
and uses the multiplication routine in the FFT program* 

This is assembled in two parts* One part is a high pass 
filter for detrending the data (removal of DC and very 
low frequencies) and the second part is for the BSD 
meas''urement* These arc assembled from 5200^^ to 5202^^ 
and 5500 ^^ to 559G^^. / 



5.2 RESULTS 


Since the accolcrometer data is not available, 
the high pass filter and integrator algorithm is tested 
for DC superimposed sine wave and the various outputs 
are given at Fig* 5*1 foi’ a section of the data. This 
program is tested for various frequencies of the input 
data and found to be satisfactory* However, the PSD 
measurement through FFT is tested using a typical track 
profile data, supplied by PDSO (Research, Development 
and Systems Organization, Lucknow) and the results are 
shown in Fig* 5*2 . 'Biis is also tested using standard 
FFT Fortran subroutines on DEC- 1 090 computer and the 
results, are found to be satisfactory* Due to the limited 
precision of signal processor which uses the 6800 micro- 
processor, the accuracy is somev/hat affected. The run 
time of double integration algorithm is approximately 
2 seconds, .v/hereas for FFT subroutine and PSD routine 
the time is approximately 18 seconds. This is mainly ■ 
due to the software multiply routine which tahes an 
approximate time of 1*67 m sec for each 16 x 16 BIT 
multiplication. A hardware multiplier using one of the 
TRW multiplier chips reduces the multiplication time by; 
a,pproxlmately 1/5} since the worst case timing including 
all the instructions and returns from the subroutine will 
be about 250 |Jl sec* 



5.3 CONCLUSIONS. 

In this thesis, a complete hardware system tor 
routine track monitoring is described in detail. All 
the individual hardware modules are tested and packed 
in a 10” x 12” x 12” hox. All the software routines 
are tested. However, the complete system performance 
is not evaluated due to the la.ck of the appropriate data 
But it is expected that the system works satisfactorily 
in the real time environment due to its consistent per- 
formance while testing individual components. 
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n 6800 


ASSEMBLY LISTING. 



* 6800 ASSEMBLY PROGRAM 

* rOR FIRST ORDER INCRE- 

* MENTAL FILTER AND 

* SIMPSON INTEGRATOR 

* FOLLOWING ARE WORKING 

* LOCATIONS FOR THE PROGRAM 


IFFC 




X! 

EQU 

SIFFC 

S8FC 




X2 

EOU 

$28FC 

3IFC 




Y1 

EOU 

$3 !FC 

A0FC 




Y2 

EQU 

$4eFC 

A905 




DELY 

EOU 

$4905 

49 10 




DELY0 

EQU 

$4910 

4915 




DELYl 

EOU 

$4915 

4920 




DELY2 

EQU 

$4920 

4925 




DELI 

EQU 

$4925 

4930 




DEL2 

EQU 

$4930 

4935 




■ INTER 

EQU 

$4935 

4940 




INTERl 

EQU 

$4940 

4945 




INTER2 

EQU 

$4945 

49 50 




TEMP 

EQU 

$4950 

4955 




TEMPI 

EQU 

$4955 

4960 




TEMP2 

EOU 

$4960 

4975 




SHICOU 

EQU 

$4975 

4980 




COUNT! 

EQU 

$4980 

4985 




C0UNT2 

EOU 

$4985 

4990 




OUT0 

EQU 

$4990 

499 5 




OUTl 

EQU 

$4995 

49A0 




0UT2 

EQU 

$49A0 





* 

MAIN 

PROGRAM START 

5000 




' 

ORG 

$5000 

5000 

7F 

IF 

FC 


CLR 

X 1 

5003 

7F 

IF 

FD 


GLR 

X 1 + 1 

5006 

7F 

IF 

FE 


CLR 

XI +2 

5009 

7F 

IF 

FF 


CLR 

XI +3 

500C 

7F 

28 

FC 


CLR 

X2 

S00F 

7F 

28 

FD 


CLR 

X2+1 

5012 

7F 

28 

FE 


CLR 

X2+2 

50 1 5 

7F 

28 

FF 


CLR 

X2 + 3 

5018 

7F 

31 

FC 


GLR 

Y1 

50 IB 

7F 

31 

FD 


CLR 

Yl + i 

50 IE 

7F 

31 

FE 


CLR 

Yi+2 

5021 

7r 

31 

FF 


CLR 

Yl+3 

5024. 

7F 

40 

FG 


CLR 

' Y2 ': ^ 

5027 

7F 

40 

FD 


CLR 

Y2 + i 

S02A 

7F 

40 

FE 


CLR 

Y2+2: 

5020 

7r 

40 

FF 


CLR 

■ , y2 + 3. ■ 

50 30 

•jf jr 

49 

05 


CLR 

DELY 

5033 

7F 

49 

06 


CLR 

DELY+1 

5036 

7F 

49 

07 


CLR 

DELY -*-2 

5039 

7F 

49 

08 


CLR 

DELY +3 

503C 

7f 

49 

10 


CLR 

DELYB 

503F 

7f 

49 

1 1 


CLR 

DELY0+1 ' 

5042 

7F 

49 

12 


CLR 

DELY0+2 .. 

5045 

7F 

49 

1 3 


CLR 

DELY 0+3 

5048 

7F 

49 

15 


CLR 

DELYl 



S04B 

7F 

49 

16 

504E 

7F 

49 

17 

5051 

7F 

49 

18 

5054 

7F 

49 

20 

5057 

7F 

49 

21 

50 5A 

7F 

49 

22 

505D 

7F 

49 

23 

5060 

7F 

49 

25 

5063 

7F 

49 

26 

5066 

7F 

49 

27 

5069 

7F 

49 

28 

506C 

7F 

49 

30 

50 6 F 

7F 

49 

31 

5072 

7F 

49 

32 

5075 

7F 

49 

33 

50 78 

7F 

49 

35 

50 7B 

7F 

49 

' 36 

507E 

7F 

49 

37 

5081 

7F 

49 

' 38 

5084 

7F 

49 

40 

508 7 

7F 

49 

41 

508A 

7F 

49 

42 

508D 

7F 

49 

43 

5090 ' 

7F 

49 

45 

5093 

7F 

49 

46 

5096 

7F 

49 

47 

5099 

IF 

49 

48 

509C 

7F 

49 

90 

609F 

7F 

49 

91 


S2A2 7F 49 92 
50AS 7F 49 93 
50A8 7F 49 95 
50AB 7F 49 96 
50AE 7F 49 97 
50B1 7F 49 98 
50B4 7F 49 A0 ■ 

50B7 7F 49 A1 
50BA 7F 49 AB 
50BD 7F 49 A3 
50C0 7F 49 50 
50C3 7F 49 55 
S0C6 7F 49 60 
50C9 7F 49 75 
50CC CE IF FF OUTER 
50CF FF 49 80 
50B2 CE 32 03 r 
50D5 FF 49 85' 

. *■ - . 

50B8 FE 49 80 MAIM 
50DB 8C JF FF • 

50DE 27. 12 
50 E0 SC 28 FF' 

50E3 27 0D 

50E5 8C 31 - 

50'£8 27 6B " 

5f£A SC . 40 FFv: ' 


CLR 

DELYl +1 

CLR 

DELYl +2 

CLR 

BELYI +3 

CLR 

DELY2 

CLR 

de:ly2+i 

CLR 

BELY2+2 

CLR 

DELY2+3 

CLR 

DELI 

CLR 

DELI +1 

CLR 

DEL 1+2 

CLR 

DELI +3 

CLR 

DEL2 

CLR 

DEL2+1 

CLR 

DEL2+2 

CLR 

DEL2+3 

CLR 

INTER 


CLR IMTER+1 

CLR ■ IMTER+2 


CER INTER+3 

CLR INTER 1 

CLR lMTERl+1 

CLR INTER 1+2 

CLR IMTERl+3 

CLR INTER2 

CLR INTER2+1 

CLR INTER2+2 

CLR IMTER2+3 

CLR OUT0 

CLR OUT0+1 

CLR OUT0+2 

CLR OUT0+3 

CLR OUTl 

CLR OUTl+1 

CLR GUT I +2 

CLR OUTl +3 

CLR 0UT2 

CLR 0UT2+1 


CLR 

0UT2+2 

CLR 

0UT2+3 

CLR 

TEMP 

CLR 

TEMP] 

CLR 

TEMP 2 

CLR 

SH I CO U 

LDK 

#£1FFF 

STX . . . 

COUNT! 

LDX 

r$3203 

STX'' ^ 

COUNT 2 


MAIN LOOP OF FILTER PROGRAM 


LDX 

COUNTl 

CPX 

#$I FFF 

'BEQ 

INNER! 

CPX ■ 

ilOSFF 

8EQ■■ 

INNERI 

CPX - 

• '"• #S'3:iFF' 

BEQ 

, INNERP 

CPX 




50 ED 

27 

66 



BE8 


INMEH2 

50 EF 

7E 

51 

B5 


JMP 


INNER 





J$5 

TAKING 

CARE OF DELAYED VALUES 






Fl3 THE 

FIRST LOOP 

50F2 

B6 

49 

43 

imtRi 

LDA 

A 

INTER 1+3 

S0F5 

B7 

49 

38 


STA 

A 

lNTER+3 

. 50F8 

E6 

49 

42 


LDA 

A 

lNtERI+2 

50F8 

B7 

49 

37 


STA 

A 

lNTER+8 

50FE 

B6 

49 

41 


LDA 

A 

INTER! +1 

5101 

E7 

49 

36 


STA 

A 

lNTER+1 

5104J 

B6 

49 

40 


LDA 

A 

INTER! 

5107 

B7 

49 

35 


STA 

A 

INTER 

510A 

B6 

49 

18 


LDA 

A 

DELYi +3 

S10D 

B7 

49 

08 


STA 

A 

DELY+3 

5110 

B6 

49 

17 


LDA 

A 

DELYI +2 

5U3 

B7 

49 

07 


STA 

A 

DELY+2 

51 16 

E6 

49 

16 


LDA 

A 

DELYI +1 

5119 

• B7 

49 

06 


STA 

A 

DELY+1 

5UC 

B6 

49 

15 


LDA 

A 

DELY 1 

51 IF 

B7 

49 

05 


STA 

A 

DELY 

5128 

B6 

49 

28 


LDA 

A 

DELI +3 

5125 

B7 

49 

13 


STA 

A 

DELY0+3 

5128 

B6 

49 

27 


L DA 

A 

DELI +2 

512B 

B7 

49 

12 


STA 

A 

DELY0+2 

512E 

B6 

49 

26 


LDA 

A 

DELl+l 

5131 

B7 

49 

1 I 


STA 

A 

DELY0+1 

513A 

B6 

49 

25 


LDA 

A 

DELI 

5137 

B7 

49 

10 


STA 

A 

DELY0 

51 3A 

B6 

49 

98 


LDA 

A 

OUTl +3 

51 3D 

B7 

49 

93 


STA 

A 

OUT0+3 

‘5140. 

B6 

49 

97 


LDA 

A 

OUTl +2 

5143’ 

B7 

49 

92 


STA 

A 

OUT0+2 

5146 

B6 

49 

96 


LDA 

A 

OUTl+1 

5149 

B7 

49 

91 


STA 

A 

OUT0+1 

514C 

B6 

49 

95 


LDA 

A 

OUTl 

514F 

b7 

49 

90 


STA 

A 

OUT0 

5152- 

7E 

51 

B5 


JMP 


INNER 






TAKING 

CARE OF DELAYED VALUES 






FOR 

THE SECOND LOOP 

5155 

B6 

49 

48 

imm2 

LDA 

A 

INTER2+3 

5158 

B7 

49 

38 


STA 

A 

INTER+3 

515B ' 

B6 

49 

47 


LDA 

A 

INTER2+2 

51 5E 

B7 

49 

37 


STA 

A 

INTER+2 

5161 

B6: 

49 

46 


LDA 

A • 

INTER2+1 

5164 

B7 

49 

36 


STA 

a: : 

IMTER+l ■ 

5167 

B 6 

49 45 


LDA A 

INTERG 

5i6A 

B7 

49 

35 


STA 

A 

INTER 

51 6 D 

B6 

49 

23 


LDA 

A 

DELYe + 3''. 

5170 

B7 

49 

08 


STA 

A 

DELy+3 

5173 

B6 

49 

22 


LDA 

A 

■DELY2+.2 ,, 

51-76 

B7 

49 07 


STA A 

cELy+2 

5179 

B6 

49 

21 


LDA 

A 

. DELY2 + I:',.':.,: 

51 7C 

87 

49, 

06 


STA 

A 

DELY+T-V^ 

517r 

E6 

49 

20 


LDA 

■ P ' 

■:0 

■ >'|)ELY2 . 

'51 82 ■ 

B7 

49 

05 : 


STA 


: UFLY.'- 

•: 5185 

B6 

49 

33 : 


LDA 

A 




5188 Bt 49 13 
S18B B6 49 32 
518E B7 49 12 
5191 B6 49 31 
5194 B7 49 11 
5197 B6 49 30 
519A 87 49 10 . ' 

5190 B6 49 A3 
51A0 B7 49 93 
51A3 B6 49 A2 
S1A6 B7 49 92 
51A9 86 49 A1 

51AC B7 49 91 
51AF B6 49 A0 
5182 B7 49 90 

5185 B6 49 38 INNER 

5 IBS F6 49 37 

51BE B0 49 08 

SIBE F2 49 07 

51.C1 B7 49 38 

51C4 B7 49 53 

51C7 F7 49 37 

51CA F7 49 52 

51CD B6 49 36 

5100 F6 49 35 

5103 B2 49 06 

5106 F2 49 05 

5109 B7 49 36 

51 DC B7 49 51 

5IDF F7 49 35 

51E2 FT 49 50 

51E5 86 05 ' 

51 E7 B7 49 75 
51EA 77 49 50 LOOP 
51E0 76 49 51 

51F8 76 49 52 

51F3 76 49 53 

51 F6 7A 49 75 
51F9 70 49 75 

51FC 26 EC 
5‘lFE A6 04 
5200 E6 03. 

5202 A0 00 

5204 09' 

5205 E? 00 

5207 B7 49 08 
520A F7 49 07 
5200 09 

520E A6 04 
5210 E6 m 
52 12 A2 00 ■ 

5214 09 . : 

5215 E2 00 

5217 . 08' - ^ 

5218 08 


STA A DELY0+3 

LDA A. m:L2*2 

StA A DELY0+2 

LDA A 0EL2+1 

STA A DELY0+1 

LDA A 0EL2 

STA A DELYe 

LDA A 0UT2+3 

STA P. OUTe + 3 

LDA A 0UT2+2 

STA A OUT0+2 

LDA A 0UT2+1 

STA A OUT0+1 

LDA A 0UT2 

STA A OUT0 

INNER LOOP OF THE PROGRAM 

LDA A INTER+3 

LDA B INTER+2 

SUE A DELY+3 

SBC B DELY+2 

STA A INTER+3 

STA A ■ TEMP+3 

STA B INTER+2 

STA B TEMP+2 

LDA A INTER+1 

LDA B INTER 

SBC A DELY+l 

SBC B DELY 

STA A iNTER+l 

STA A TEMP+1 

STA B INTER 

STA B . TEMP 

LDA A #S05 

STA A SHICOU 

ASR TEMP 

ROR TEMP+1 

ROR TEMP+8 

ROR TEMP+3 ■ 

DEC SHICOU 

TST SHICOU 

BNE LOOP 

LDA A 4^X • 

LDA B 3^X 


SUB 

DEX 

A 

■ 0aX , : : 

SBC 

B 

, 

STA 

A 

. : DELY + 3 

STA 

LEX 

B 

DELY +2 

LDA 

A 

4#X 

LDA 

B 

3#X 

SBC 

LEX 

A 

0»X 

SBC 

I NX 
INX 

B 

0#X 



5SI9 

08 



INK 



5SIA 

B7 

49 

06 

STA 

A 

DELY-f-l 

5210 

F7 

49 

05 

STA 

B 

DELY 

5220 

B6 

49 

08 

LDA 

A 

DELY<-3 

5223 

F6 

49 

07 

L0A 

B 

DELY +2 

5226 

BB 

49 

53 

ADC 

A 

TEMPOS 

5229 

F9 

49 

52 

ADC 

B 

TEMP+2 

522G 

B7 

49 

08 

ETA 

A 

DELY+3 

522F 

F7 

49 

07 

StA 

B 

DELY -*-2 

5232 

B6 

49 

06 

LDA 

A 

D£LY+i 

5235 

F6 

49 

05 

LDA 

B 

DELY 

5238 

E9 

49 

51 

ADC 

A 

TEMP+I 

523B 

F9 

49 

50 

ADC 

B 

TEMP 

523E . 

B7 

49 

06 

STA 

A 

DELY 4- S 

5241 

F7 

49 

05 

STA 

B 

DELY 

5244 

B6 

49 

08 

LDA 

A 

DELY+3 

5247 

F6 

49 

07 

LDA 

B 

DELY+2 

524A 

BB 

49 

13 

ADD 

A 

DELY0+3 

52 4 D 

F9 

49 

12 

ADC 

B 

DELY0+2 

5250 

B7 

49 

13 

STA 

A 

DELY0+3 

5253 

B7 

49 

58 

STA 

A 

TEMP 1 + 3 

5256 

F7 

49 

12 

STA 

B 

DELY0+2 

5259 

F7 

49 

57 

STA 

B 

TEMP 1+2 

525C 

B6 

49 

06 

LDA 

A 

DELY+l 

525F 

F6 

49 

05 

LDA 

B 

DELY 

5262 

B9 

49 

1 1 

ADC 

A 

DELYB+l 

5265 

F9 

49 

10 

ADC 

B 

DELY0 

5268 

B7 

49 

11 

STA 

A 

DELY0+1 

526B 

B7 

49 

56 

STA 

A 

TEMF1+ 1 

526E 

F7 

49 

10 

STA 

B 

DELY0 

5271 

F7 

49 

55 

STA 

B 

TEMPI 

5274 

86 

04 


LDA 

A 

#$04 

5276 

B7 

49 

75 

STA 

A 

SHICOU 

5279 

77 

49 

55 LOOPl 

ASR 


TEMPI 

527C 

76 

49 

56 

ROE 


TEMPI + l 

52 7F 

76 

49 

57 

ROR 


TEMP 1+2 

5282 

76 

49 

58 

ROE 


TEMPI +3 

5285 

7A 

49 

75 

DEC 


SHICOU 

5288 

70 

49 

75 

TST 


SHICOU 

528B 

26 

EC 


BNE 


LOOP! 

5280 

fjr 

49 

60 

STX 


TEMP2 


* BEANCHING TO INTEGRATION EOUTIME 


5290 

BB 

'54 

46 

JSR 

INTEGR 

5293 ■ 

■ FE 

49 

85 

■LBX'.''.' 

COUNTS 

5296 

B6 

49 

58 

. .. LBA A . 

TEMP 1 + 3 

5299 

A7 

m 


STA A 

0>X 

.52 9B 

09 



. ■ DEX 


529C 

B6 

49 

57 

LDA A 

TEMPI +2 

529F- 

A7 

00 


STA A 

0#X 

52A1 

09 



DEX 


52A2 

B& 

49 

56 

LDA A 

TEMPi + 1 

S2A5 

A7 

00 


STA A 

0^X. 

52A7 

09 



DEX 


52A8 

B6 

49 

55“' 

LDA A 

■ TEMPI-:: 

52AB 

A7 

00 


STA A 

0#X 

S2AD 

FE 

49 

60 

LDX 

TEMP2: 



52B0 . 

■ EC 

S7 

FB 


CPX 

#S27FB 

52B.3 

27 

1 2 



e£@ 

STOFEl 

5SB5 

SC 

3i 

FB 


CPX 


52B8 

2? 

23 



BEQ 

ST0R2 

5 2 BA 

8C 

39 

FB 


CPX 

#$39FB 

S28B 

27 

2i 



BEQ 

StORS 

S2BF 

8C 

48 

FB 


CPX 

#$4aFB 

S2CS 

27 

ID 



BEE 

ST0R6 

52C4 

7E 

S4 

00 


JMP 

NEXT 






TAKSKG 

CARE OF LAST 





* 

LOCATIO^^ VALUES 

52C7 

A6 

04 


StOHES 

LDA A 

4# X 

52G^ 

37 

28 

ff' 


STA A 

X2 + 3 

52CG 

A6 

03 



/a 

3^X 

52CE 

B7 

28 

F£ 


STA A 

X2->2 

S2DI 

A6 

02 



L DA A 

2^X 

52D3 

B7 

28 

FD 


STA A 

X2-9-1 

52D6 

A6 

01 



LDA A 


52D8 

37 

28 

FC 


STA A 

X2 

52DB 

20 

06 



BRA 

PRSTOl 

52DD 

20 

71 


ST0R2 

BRA 

ST0RE2 

52DF 

20 

&A 


ST0R5 

BRA 

STORES 

52E1 

20 

6A 


ST0R6 

BRA 

STORES 

52E3 

B6 

49 

08 

PRSTOl 

LDA A 

DELY+3 

■52E6 

B7 

49 

» 18 


STA A 

DELYI-s-a 

52E,9 

B6 

49 

' 07 


LDA A 

DELY+2 

52EC 

B7 

49 

17 


STA A 

DELYl+2 

52EF 

B6 

49 

06 


LDA A 

DELY+! 

52F2 

B7 

49 

16 


STA A 

DELYl+1 

52F5 

B6 

49 

05 


LDA A 

DELY 

52F8 

B7 

49 

15 


STA A 

DELY 1 

52FB 

B6 

49 

38 


LDA ' A 

INTER+3 

52FE 

B7 

49 

43 


STA A 

INTER! +3 

5301 

B6 

49 

37 


LDA A 

INTER+2 

5304 

B7 

49 

42 


STA A 

INTER] +2 

5307 

B6 

49 

36 


LDA A 

XNTER+ 1 

530A 

B7 

49 

4! 


STA A 

INTER! +1 

530D 

B6 

49 

' 35 


LDA A 

INTER 

5310 

B7 

49 

40 


STA A 

INTER! 

5313 

B6 

49 

13 


LDA A 

DELY0+3 

5316 

B7 

49 

28 


STA A 

DELI +3 

5319 

B6 

49 

' 12 


LDA A 

DELY0+2 

531C 

B7 

49 

27 


STA A 

DELI +2 ■ 

53 IF 

B6 

49 

il 


LDA A 

DELY 0+1 

5322 

B 7 

49 

' 26 


STA A 

DELl+l 

5325 

B6 

49 

10 


LDA A 

DELY0 

5328 

B7 

49 

25 


STA A 

DELI 

532B 

B6 

49 

93 


LDA A 

OTJT0+3 ' 

532E 

F6 

49 

92 


lda‘ B 

•OUT0+2 

5331 

B7 

49 

98 


STA A 

OUTl+3 

5334 

F7 

*49 

97 


STA B 

OUTl+2 

5337 

B6 

49 

91 


LDA A 

OUT0+1 

533A 

F6 

49 

90 


LDA B 

OUT0 

S33D 

B7 

49 

96 


STA A 

OUTl+l 

5340 

F7 

49 

95 


STA B 

OUTl 

5343 

ar 


ft , 

, 

m 




5346 

2? 

iF 



BE© 


FOLLOI 

5348 

7E 

54 

17 


JMP 


S ECO ML 

534E 

20 

ID 


STORES 

BHA 


STORES 

5341; 

7E 

S3 

E6 

STORES 

JMP 


ST0RE4 

5350 

A6 

04 


STORES 

LDA 

A 

4#X 

S352 

B7 

IF 

f-f 


STA 

A 

Xl+3 

5355 

A6 

03 



LDA 

A 

3*X 

535? 

B7 

IF 

FE 


STA 

A 

XI +2 

S35A 

A6 

02 



LDA 

A 

2,X 

535C 

B7 

IF 

F,D 


■ STA 

A 

Xi + 1 

53 5F 

A6 

01 



LDA 

A 

l^x 

536 i 

S7 

IF 

FC 


STA 

A 

xi 

5364 

7E 

52, 

E3 


JMP 


PRSTO 1 

5367 

?E 

54 

26 

FOLLOI 

JMP 


THIRD 

536A 

A6 

04 


ST0HE3 

LDA 

A 

4^ X 

536C 

B? 

40 

FF- 


STA 

A 

V2 + 3 

536F 

A6 

03 



LDA 

A 

3*X 

537 J 

B? 

48 

FE 


STA 

A 


5374 

A6 

02 



LDA 

A 

2^ X 

5376 

B7 

40 

FD 


STA 

A 

V2+1 

5379 

A6 

01 



LDA 

A 

I 

537E 

B7 

40 

FC 


STA 

A 

Y2 

537E 

. B6 

49 

03 

PRST02 

LDA 

A 

DELY+3 

538 i 

B7 

49 

23 


STA 

A 

DELY2+3 

5384 ' 

B6 

49 

07 


LDA 

A 

DELY+2 

5387 • 

B7 

49 

22 


STA 

A 

DELY2^-2 

538A 

B6 

49 

06 

* 

LDA 

A 

DELY+I 

538D 

B7 

49 

21 


STA 

A 

DELY2+I 

5390 

B6 

49 

05 


LDA 

A 

DELY 

5393 

B7 

49 

20 


STA 

A 

DELY2 

5396 

B6 

49 

38 


LDA 

A . 

lNTER+3 

5399 

B7 

49 

48 


STA 

A 

1NTER2+3 

539C 

B6 

49 

37 


LDA 

A 

INTER+2 

539F 

E7 

49 

47 


STA 

A 

INTER2+2 

53A2 

B6 

49 

36 


LDA 

A 

INTER+1 

53A5 

B7 

49 

46 


STA 

A 

INTER2+1 

53A8 

B6 

49 

35 


LDA 

A 

INTER 

S3AB 

B7 

49 

45 


STA 

A 

INTER2 

53At 

B6 

49 

13 


LDA 

A 

DELY0+3 

53BI 

B7 

49 

33 


STA 

A 

DEL2+3 

53B4 

B6 

49 

12 


LDA 

A " . 

DELY0+2 

53B7 

B7 

49 

32 


STA 

A 

CEL2+2 

S3BA 

B6 

49 

I 1 


LDA 

A 

DELY0+I 

5300 

D7 

49 

31 


STA 

A 

DEL2-H 

5 3 C 0 

B6 

49 

10 


LDA 

A 

DELY0 

53C3 

.07 

49 

30 . 


STA 

A 

DEL© . 

53C6 

B6 

49 

.93 


LDA 

A 

OUT0+3 

53C9 

B7 

49 

A3 


STA 

A 

0UT2f3 

S3CC 

B6 

49 

92 


LDA 

A 

0UT@*2 

S3CF 

B? 

49 

A2 


STA 

A 

0UT2+2 

53D2 

86 

49 

91 


LDA 

A 

OUT0+I 

53D5 

B7 

49 

At 


STA 

A 


S3C8 

&6 

.49 

90 


LI A 

A 

om# 

53D0 

B7 

49 

M 


STA 

A 

OUT 2 

sat't 

SC 

48 

FB 


CPK 


# i43Fb 

S3fci 

27 

lA 



SEQ 


. rOLLO? 







^ 3 E 3 

'?E 

54 

33 

5.3E6 

A6 

i4 


S3E8 

B? 

31 

Ff 

53ES 

A6 

03 


S3ED 

ST 

31 

TE 

53F0 

A6 

02 


33F2 

37 

3 1 

FD 

S3FS 

Ae 

Bl 


: 33F7 

S'f 

3 1 

FC 

S3FA 

7£ 

53 

r & 

53FD 

7E 

54 

44 

S 4 0 

FE 

49 

80 

3 4 2 3 

38 



5404 

08 



34 f 5 

S8 



S4 0 5 

08 



34g? 

FF 

49 

BB 

540 A 

FS. 

4 9. 

BS 

S4fiD 

§8 : 



S40E 

0S 



340F 

08 



5411 

08 



S 4 J ^ 

FF 

49 

05 

54 1 4 , 

?E 

5i 

0$ 


tl p 

^ a .^1 

S T A A 

4# i*.^ lA ^»i 

ST A A 

'r< A # 

^ s% 

. STA A 
^ &} A A 
STA A 
JM? 

AT EXT LDX 

XPJX 
I&"X 

r.«x 

■ iMX ' 

'■■ .SIX 
LBX 
ISX . 

mx 

INX 

ISX 

%J'' 

*« 


F O'UFiTH 
X 

'/I 3 . 

3 j» X 


Yt ■*■2 
2^X 
Y 2 * T 
ji X 
Y1 

FRSTC2 
GUT . 


COUMTi 
COURT 2 


C0WT8 

MAIM 








CE 3J FF S2C0KC LDX 
FF 49 80 . STX 


■^s3irF 

Tiiiilir 

illi'iliS 

■tliliii: 

^iiwiii 

•iiliiS" 

iiiaiA::-':: 

iiifiis 

caus'c2 


imm 








5^# If 
, ^ 4 ' ii A 
^S4iS 
i4:ii :::ArP 
S4i'l 



546E 

,F7 

49 

90 

STA B 

01)10 

5471 

B7 

49 

56 

STA A 

TEHPl+1 

5474 

F7 

49 

55 

STA B 

TEKP! 

5477 

39 



RTS 


5478 




END 



STATEME;MTS *453 
FREE BYTES » 19460 
NO ERRORS DETECTED 



i 





* ' 

FFT IN 

PLACE PROGRAM 




* 

N IS NO 

•OF SAMPLES 




* 

N SHOULD BE POWER OF TWO 

0080 




ORG 

$80 




* 

WORKING LOCATIONS FOR 




* 

MULTIPLICATION 




* 

USING BOOTHS ALGORITHM 

0080 



Y 

RMB 

2 

0082 



XX 

RMB 

2 

0084 



u 

RMB 

4 

0088 



FF 

RMB 

1 

4000 




ORG 

$4000 

4000 

01 


FFTIN 

NOP 


4001 

7E 

43 74 


JMP 

FFT 




* 

257 VALUES OF SINE ARE 




* 

STORED 

FROM 0 TO PIE/2 




* 

INSTEPS 

OF 2PIE/ 1024 

4004 

00 

00 

SIN 

FDB 

$0000 

4006 

00 

C9 


FDB 

$00C9 

4008 

01 

92 


FDB 

$0192 

400A 

02 

5B 


FDB 

S025B 

400C 

03 

24 


FDB 

$0324 

40 0E 

03 

FD 


FDB 

$03FD 

4010 

04 

B6* 


FDB 

$04B6 

4012 

05 

7F. 


FDB 

$057F 

4014 

06 

48 


FDB 

$0648 

4016 

07 

1 1 


FDB 

$0711 

4018 

07 

D9 


FDB 

$07D9 

401 A 

08 

A2 


FDB 

$08A2 

40 1C 

09 

6 B 


FDB 

$096B 

401E 

0A 

33 


FDB 

S0A33 

4020 

0A 

FB 


FDB 

$0AFB 

4022' 

0B 

C4 


FDB 

$0BC4 

4024 

0C 

8C 


FDB 

$0G8C 

4026 

0D 

54 


FDB 

$0D54 

4028 

0E 

1C 


FDB 

S0E1C 

402A 

0E 

E4 


FPB 

$0EE4 

402C 

0F 

AB 


' FPS ' . 

' $0FAB 






404E 

IC 

D0 

FDB 

$1CD0 

4050 

ID 

93 

FDB 

S1D93 

4052 

IE 

57 

FDB 

$1E57 

4054 

IF 

lA . 

FDB 

SIFIA 

4056 

IF 

DD 

FDB 

SIFDD 

4058 

20 

9F 

FDB 

S209F 

405ft 

21 

62 

FDB 

S2162 

405C 

22 

24 

FDB 

$2224 

40 5 E 

22 

E5 

FDB 

$22E5 

4060 

23 

A7 

FDB 

$23A7 

4062 

24 

67 

FDB 

$2467 

4064 

25 

28 

FDB 

$2528 

4066 

25 

E8 

FDB 

$25E8 

4068 

26 

A8 

FDB 

$26A8 

406ft 

27 

68 

FDB 

$2768 

40 6 C 

28 

27 

FDB 

$2827 

406E 

28 

E5 ' 

FDB 

$28E5 

4070 

29 

A4 

FDB 

$29A4 

4072 

2A 

62 

FDB 

$2A62 

4074 

2B 

IF 

FDB 

$2B1F 

.4076 

2B 

DC 

FDB 

$2BBC 

4078 

2C 

99 

FDB 

$2C99 

407ft 

2D 

55 

FDB 

$2D55 

407C 

2E 

1 1 

FDB 

$2E1 1 

407E 

2E 

CC 

FDB 

$2ECC 

4080 

2F 

87 

FDB 

$2F87 

4082 

30 

42 

FDB 

$3042 

4084 

30 

FC 

FDB 

$30FC 

4086 

31 

B5 

FDB 

$31B5 

4088 

32 

6E 

FDB 

$326E 

408ft 

33 

27 

FDB 

$3327 

408 C 

33 

DF 

FDB 

$33DF 

408E 

34 

9 7 

FDB 

$3497 

4090 

35 

4E 

FDB 

$354E 

4092 

36 

04 

FDB 

$3604 

4094 

36 

BA 

FDB 

$36 BA 

4096 

37 

70 

FDB 

$3770 

4098 

38 

25 

FDB 

$3825 

409A 

38 

D9 

FDB 

$38 D9 

409C 

39 

8D 

FDB 

$398D 

40 9E 

3A 

40 

FDB 

$3A40 

40 A0 

3A 

F3 

FDB 

$3AF3 

40A2 

3B 

A5 

FDB 

$3BA5 

40A4 

3C 

57 

FDB 

$3C57 

40 A6 

3D 

08 

FDB 

$3D08 

40A8 

3D 

B8 

FDB 

$3DB8 

40 Aft 

3E 

68 

FDB 

$3E68 

40AC 

3F 

1 7 

FDB 

$3F17 

40AE 

3F 

C6 

FDB 

$3FC6 

40B0 

40 

74 

FDB 

$4074 

40B2 

41 

21 

FDB 

$4121 

40B4 

41 

CE 

FDB 

$4lCE 

40B6 ’ 

42 

7ft 


; 

40B8 

43 

26 

filly; 


40BA 

43 

D1 • 



40BC 


7B 





40BE 

45 

24 

FDB 

S4524 

40 C0 

45 

CD 

FDB 

$45CD 

40C2 

46 

75 

FDB 

S4675 

40C4 

47 

ID 

FDB 

$471D 

40C6 

47 

C4 

FDB 

$47C4 

40C8 

48 

6A 

FDB 

$486A 

40CA 

49 

0F 

FDB 

$490F 

40CC 

49 

B4 

FDB 

$49B4 

40CE 

4A 

58 

FDB 

S4A58 

40D0 

4A 

F -B 

FDB 

$4AFB 

40 D2 

4E 

9E 

FDB 

$4E9E 

40D4 

4C 

40 

FDB 

S4C40 

40D6 

4C 

El 

FDB 

♦ S4CE1 

40D8 

4D 

81 

FDB 

S4D81 

40 DA 

4E 

21 

FDB 

J4E21 

40DC 

4E 

C0 

FDB 

S4EC0 

40DE 

4F 

5E 

FDB 

$4F5E 

40E0 

4F 

FB 

FDB 

$4FFB 

40E2 

50 

98 

FDB 

$5098 

40E4 

51 

34 

FDB 

$5134 

40E6 

51 

CF 

FDB 

$51CF 

40 E8 

52 

69 

FDB 

$5269 

40 EA 

53 

03 

FDB 

$5303 

40EC 

53 

9B 

FDB 

$539B 

40EE 

54 

33 

FDB 

$5433 

40F0 

54 

CA 

FDB 

$54CA 

40F2 

55 

60 

FDB 

$5560 

40F4 

55 

F6 

FDB 

$55F6 

40F6 

56 

8A 

FDB 

. $56 8 A 

40F8 

57 

IE 

FDB 

$57 IE 

40FA 

57 

B1 

FDB 

$57B1 

40FC 

58 

43 

FDB 

$5843 

40FE 

58 

D4 

FDB 

$58D4 

4100 

59 

64 

FDB 

$5964 

4102 

59 

F4 

FDB 

$59F4 

4104 

5A 

82 

FDB 

$5A82 

4106 

5B 

10 

FDB 

$5B10 

4108 

5B 

9D 

FDB 

$5B9B 

410A 

5C 

29 

FDB 

$5C29 

4 1 0 C 

5G 

B4 

FDB 

$5CB4 

4 1 0E 

5D 

3E 

FDB 

$5D3E 

4110 

5D 

C8 

FDB 

$5 DC 8 

4112 

5E 

50 

FDB 

$5E50 

41 14 

5E 

D7 


$5 ED 7 

4116 

5F 

5E 

FDB 

$5F5E 

4118 

5F 

E4 

FDB 

$5FE4 

41 1 A 

60 

68 

FDB 

$606 8 

41 1C 

60 

EC 

FDB 

$60EC 

41 IE 

61 

6F 

FDB 

$616F 

4120 

61 

FI 

FDB 

$61F1 

4122 

62 

72 

FDB 

$6272 

4124 

62 

F2 

FDB 

$62F2 

4126 

63 

7 1 

, FDB 


4128 

63 

EF 



412A 

64 

6C 


$646C 

412C 

64 

E9 

FDB 

$64E9' 



412E 

65 

6 4 

FDB 

$6564 

4130 

65 

DE 

FDB 

$65DE 

4132 

66 

57 

FDB 

$6657 

4134 

66 

D0 

FDB 

$66D0 

4136 

67 

47 

FDB 

$6747 

4138 

67 

BD 

FDB 

$67BD 

413A 

68 

32 

FDB 

$6832 

413C 

68 

A7 

FDB 

$68A7 

413E 

69 

lA 

FDB 

$691A 

4140 

69 

8C 

FDB 

$698C 

4142 

69 

FD 

FDB 

$69FD 

4144 

6A 

6E 

FDB 

$6A6E 

4146 

6A 

DD 

FDB 

$6ADD 

4148 

6B 

4B 

FDB 

$6B4B 

414A 

6B 

B8 

FDB 

$6BB8 

414C 

6C 

24 

FDB 

$6C24 

414E 

6C 

8F 

FDB 

$6C8F 

4150 

6C 

F9 

FDB 

$6CF9 

4152 

6D 

62 

FDB 

$6D62 

4154 

6D 

CA 

FDB 

$6 DCA 

4156 

6E 

31 

FDB 

$6E31 

4158 

6E 

97 

FDB 

$6E97 

415A 

6E 

FB 

FDB 

$6EFB 

41 5C 

6F 

5F 

FDB 

$6F5F 

415E 

6F 

C2 

FDB 

$6FC2 

4160 

70 

23 

FDB 

$7023 

4162 

70 

83 

FDB 

$7083 

4164 

70 

E3 ’ 

FDB 

$70E3 

4166 

71 

41 

FDB 

$7141 

4168 

71 

9E 

FDB 

$719E 

416A 

71 

FA 

FJDB 

$71 FA 

416C 

72 

55 . 

FDB 

$7255 

416E 

72 

AF 

FDB 

$72AF 

4170 

73 

08 

FDB 

$7308 

4172 

73 

5F 

■ 'I'i^ 


1174 

73 

B6 

FDB 

$73B6 

4176 

74 

0B 

FDB 

$740B 

4178 

74 

60 

FDB 

$7460 

417A 

74 

B3 

FDB 

$74B3 

417C 

75 

05 

FDB , 

$7505 

417E 

75 

56 

FDB , 

$7556 

4180 

75 

A6 

FDB 

$75A6 

4182 

75 

F 4 

■ vFDB ■ 

$75F4 

4184 

76 

42 

FDB 

$7642 

4186 

76 

8E 

FDB 

$768E 

4188 

76 

D9 

FDB 

$76D9 

418A 

77 

23 

FDB 

$7723 

418C 

77 

6C 

FDB 

$776C 

41 8E 

77 

B4 

FDB 

$77B4 

4190 

77 

FB 

FDB 

$77FB 

4192 

78 

40 

FDB 

$7840 

4194 

78 

85 

FDB 

$7885 

4196 

78 

C8 

FDB 

$78C8 

4198 

79 

0A 

FDB 


41 9A 

79 

4A 

FDB 


419C 

79 

8A 

FDB ^ 




420A 

00 

00 

N 

FDB 

$0000 




♦ 

lA INITIAL VALUE IS N/2 

420C 

00 

00 

lA 

FDB 

$0000 




* 

LL STORES L0G2N 

420E 

00 

00 

LL 

FDB 

$0000 

4210 

00 

00 

HH 

FDB 

$0000 ; 




* 

FOLLOWING PARAMETERS SET 




* 

OFFSETS FOR CALCULATING 




* 

SINE j 

AND COSINE VALUES 

4212 

02 

00 

07512 

FDB 

$0200 

4214 

00 

00 : 

S1024 

FDB 

$0000 

4216 

04 

00 

C1024 

FDB 

$0400 

4218 

Wit 

0k3 

SI512 

FDB 

$0000 

42 lA 

00 

00 

SN512 

FDB 

$0000 

42 1C 

FE 

00 

C512 

FDB 

$FE00 ' 




♦ 

INNER 

LOOP PARAMETERS 

42 IE 

00 

01 

IB 

FDB 

$000 1 1 

4220 

00 

00 

IC 

FDB 

$0000 1 

4222 

00 

00 

IB 

FDB 

$0000 ; 

4224 

00 

00 

I2D 

FDB 

$0000 ; 

4226 

00 

00 

K 

FDB 

$0000 1 




* 

SINE 

IS TEMPORARY 




* 

LOCATION TO SINE 

4228 

00 

00 

SINE 

FDB 

$0000 ; 




♦ 

CO SIN 

IS TEMPORARY 




*■ 

LOCATION TO COSINE 

422A 

00 

00 

COSIN 

FDB 

$0000 ; 

422C 

00 


COUNT 

FCB 

$00 




♦ 

FLAGS 

FOR TESTING 

422D 

00 


FLAGl 

FCB 

$00 1 

422E 

00 

1 

FLAG 2 

FCB 

$00 

422F 

00 


FLAG 

FCB 

$00 

4230 

00 

00 

EL 

FDB 

$0000 

42 32 

00 

00 

EH 

FDB 

$0000 

42 34 

00 

00 

DL 

FDB 

$0000 

42 36 

00 

00 

DH 

FDB 

$0000 




* 

REAL 

PART CURRENT 




* 

ADDRESS STORE 

4238 

00 

00 

MIR 

FDB 

$0000 ; 




* 

I M AG . 

PART CURENT 




♦ 

- ADDRESS STORE 

42 3A 

00 

00 

Mil 

FDB 

$0000 1 

42 3 C 

00 

00 

NMR 

FDB 

$0000 

42 3 E 

00 

00 

NMI 

FDB 

$0000 




^ * 

CURRENT REAL PART STORE 

4240 

00 

00 

XYRMI 

FDB 

$0000 1 




* 

CURRENT IMAG.PART STORE 

4242 

00 

00 

XYIMI 

FDB 

$0000 } 




* 

WORKING LOCATIONS TO 




♦ 

STORE 

REAL* IM AG. PARTS 

4244 

00 

00 

AREAL 

FDB 

$0000 ; 

4246 

00 

00 

AIMAG 

FDB,- 

$0000 ; 

4248 

00 

00 

REALl 

; FDR 

$0000 

42 4A 

00 

00 

BREAL 



424C 

00 

00 

IMAGl 



42 4E 

00 

00 

BIMAG 


$0000 


420A 

00 

00 

N 

FDB 

$0000 




+ 

I A INITIAL VALUE IS N/2 

420C 

00 

00 

lA 

FDB 

$0000 




* 

LL STORES L0G2N 

420E 

00 

00 

LL 

FDB 

$0000 

4210 

00 

00 

HH 

FDB 

$0000 ; 




+ 

FOLLOWING PARAMETERS SET 




* 

OFFSETS FOR CALCULATING 




♦ 

SINE AND COSINE VALUES 

4212 

02 

00 

CT512 

FDB 

$0200 

4214 

00 

00 

SI 024 

FDB 

$0000 

4216 

04 

00 

C1024 

FDB 

$0400 

4218 


0I<3 

SI512 

FDB 

$0000 

421A 

00 

00 

SN51 2 

FDB 

$0000 

42 1C 

FE 

00 

C5I2 

FDB 

$FE00 ' 




* 

INNER 

LOOP PARAMETERS 

42 IE 

00 

01 

IB 

FDB 

$0001 J 

4220 

00 

00 

I C 

FDB 

$0000 1 

4222 

00 

00 

ID 

FDB 

$0000 ; 

4224 

00 

00 

I2D 

FDB 

$0000 ; 

4226 

00 

00 

K 

FDB 

$0000 ; 




* 

SINE 

IS TEMPORARY 




* 

LOCATION TO SINE 

4228 

00 

00 

SINE 

FDB 

$0000 ; 




* 

COSIN 

IS TEMPORARY 




♦ 

LOCATION TO COSINE 

422A 

00 

00 

COSIN 

FDB 

$0000 ; 

422C 

00 


COUNT 

FCB 

$00 




+ 

FLAGS 

FOR TESTING 

422D 

00 


FLAGl 

FCB 

$00 

422E 

00 

1 

FLAG2 

FCB 

$00 

422F 

00 


FLAG 

FCB 

$00 

4230 

00 

00 

EL 

FDB 

$0000 

42 32 

00 

00 

EH 

FDB ■ 

$0000 

42 34 

00 

00 

DL 

FDB 

$0000 

42 36 

00 

00 

DH 

FDB 

$0000 




■ . * ■ 

REAL 

PART CURRENT 




* 

ADDRESS STORE 

4238 

00 

00 

MIR 

FDB 

$0000 } 




* 

IMAG.: 

PART CURENT 




* 

" ADDRESS STORE 

42 3A 

00 

00 

Mil 

FDB 

$0000 ; 

42 3C 

00 

00 

NMR 

FDB 

$0000 

42 3 E 

00 

00 

NMI 

FDB 

$0000 




* 

CURRENT REAL PART STORE 

4240 

00 

00 

XYRMI 

FDB 

$0000 1 




* 

CURRENT IMAG. PART STORE 

4242 

00 

00 

XYIMI 

FDB 

$0000 ; 




*■ 

WORKING LOCATIONS TO 




* 

STORE 

REAL* IMAG. PARTS 

4244 

00 

00 

AREAL 

FDB 

$0000 i 

4246 

00 

00 

AIMAG 

FDB 

$0000 ; 

4248 

00 

00 

REALl 

, FP3 

$0 000 

42 4A 

00 

00 

BREAL 

FDB . 


424C 

00 

00 

IMAGl 

FDB ■ 


aiAi' ' 

00 

00 

BIMAG 

FDB 

$0000 



4250 

00 


FLAG4 

FCB 

$00 

42 51 

00 


BFLAG 

FCB 

$00 

4252 

00 


LL8 

FCB 

$00 




♦ 

TEMP ADDRESS LOCATIONS 




+ 

OF BIT 

REVERSAL ROUTINE 

4253 

00 

00 

TREAL 

FDB 

$0000 } 

4255 

00 

00 

TIMAG 

FDB 

$0000 ; 




♦ 

BITRM STORES BIT REVERSED# 

4257 

m 

00 

BURN 

FDB 

$0000 ; 




* 

MAJOR WORK LOCATIONS 




* 

OF FFT 

PROGRAM 

42 59 

00 


H 

FCB 

$00 

42 5A 

00 


L 

FCB 

$00 

425B 

00 


D 

FCB 

$00 

425C 

00 


E 

FCB 

$00 ; 

42 5 D 

00 


B1 

FCB 

$00 ; 

425E 

00 


C 

FCB 

$00 i 




* 

DPMUL IS A SUBROUTINE TO 




* 

MULTIPY TWO’S COMPLIMENT 




* 

NO. USING BOOTHS ALGOITHM 




* 

MULTIPLIER AND MULTIPLICAND 




* 

SHOULD 

BE STORED AT DE>HL ' 




♦ 

WORK LOCATIONS AND THE 




♦ 

RESULT 

WILL BE AT HL 

42 5F 

01 


DPMUL 

NOP 


4260 

FE 

42 5B 


LDX 

D 

4263 

DF 

80 


STX 

Y 

4265 

FE 

42 59 


LDX 

H 

4268 

DF 

82 


STX 

XX 

42 6A 

CE 

00 05 


LDX 

#5 

426D 

4F 



CLR A 


42 6 E 

A7 

83 

LPl 

STA A 

U- 1 >x 

4270 

09 



DEX 


4271 

26 

F B 


BNE 

LPl 

4273 

CE 

00 10 


LDX 

#16 

4276 

96 

81 

LP2 

LDA A 

Y+l 

4278 

84 

01 


AND A 

#1 

42 7A 

16 



TAB 


42 7B 

98 

88 


EOR A 

FF 

42 7D 

27 

ID 


BED 

SHIFT 

427F 

5D 



TST B 


4280 

27 

0E 


BEQ : 

ADD 

4282 

96 

85 


LDA A 

/'■ U+1 ■ 

4284 

D6 

84 


LDA B 

^.■U'.V- 

4286 

90 

83 


SUB A 

XX+l 

4288 

D2 

82 


SBC B 

.XX 

428A 

97 

85 


STA A 

U+1 

428C 

D7 

84 


STA B 

u 

428E 

20 

0C 


BRA 

SHIFT 

4290 

96 

85 

ADD 

LDA A 

U+1 

4292 

D6 

84 


LDA B 

u 

4294 

9B 

83 


ADD A 

XX+l 

4296 

D9 

82 


ADC B 


4298 

97 

85 



ilfciiilKI; 

42‘9A 

D7 

8 4 


STA B 


429C 

7F 

00 88 

SHIFT 

. CL^ 




429F 

76 

00 

80 

42A2 

76 

00 

81 

42A5 

79 

00 

88 

42A8 

77 

00 

84 

42AB 

76 

00 

85 

42AE 

76 

00 

86 

42B1 

76 

00 

87 

42B4 

09 



42B5 

26 

BF 


42B7 

78 

00 

86 

42BA 

79 

00 

85 

42BD 

79 

00 

84 

42C0 

96 

84 


42C2 

D6 

85 


42C4 

B7 

42 

59 

42C7 

F7 

42 

5A 

42CA 

39 




* 

* 

* 

♦ 

42CB 01 COMPL 

42CC 4F 

42CD B6 42 5A 

42D0 F6 42 59 

42D3 53 

42D4 43 

42D5 8B 01 

42D7 C9 00 

42D9 B7 42 5A 

42DC F7 42 59 

42DF 39 

* 

♦ 

42E0 01 OFL 

42E1 86 00 

42E3 B7 42 50 

42E6 B6 42 59 

42E9 84 80 

42EB 27 08 

42ED 86 01 

42EF B7 42 50 

42F2 BD 42 CB 

42 F5 01 DVIDE 

42F6 0C 

42F7 76 42 59 

42FA 76 42 5A 

42FD B6 42 50 

4300 81 01 

4302 2B 03 

43,04 BD 42 CB 

4307 39 OVER 

♦ 

♦ 

'■ ' ♦ 


ROR Y 

ROR Y+1 

ROL FF 

ASR U 

ROR U+1 

ROR U+2 

ROR U+3 

DEX 

BNE LP2 

ASL U+2 

ROL U+1 

ROL U 

LDA A U 

LDA' B U+1 

STA A H 

STA B L 

RTS 

SUBROUTINE COMPL FINDS 

TWO'S COMPLIMENT OF 16 

BIT NUMBER STORED AT HL 

RESULT AT HL 

NOP • 

CLR A 

LDA A L 

LDA B H 

COM B 

COM A 

ADD A #S01 

ADC B #$00 

STA A L 

STA B H 

RTS 

SUBROUTINE OFL DIVIDES 
16 BIT NUMBER STORED 
AT HL BY 2; RESULT AT HL , 
NOP 

LDA A #$00 

STA A FLAG4 
LDA A H 

AND A #$80 

BEQ DVIDE 

LDA A #$01 

STA A FLAG4 

JSR COMPL 

NOP 

CLC ' 

ROR H 

ROR L: 

LDA A FLAG4 

CMP A #$01 

BMI OVER 

JSR COMPL 

SUBROUTINE BITR FINDS 
BIT REVERSED NO. OF-A; 



* 

* 


NO. STORED AT HL 
RESULT AT HL 


4308 

01 



BITR 

NOP 


4309 

FE 

42 

59 


LDX 

H 

430C 

FF 

42 

5B 


STX 

D 

430F 

B6 

4£ 

! 51 


LDA A 

BFLAG 

4312 

81 

01 



CMP A. 

#$01 

4314 

2B 

3A 



BMI 

LLLE8 

4316 

B6 

42 

5C 

LLGT8 

LDA A 

E 

4319 

C6 

08 



LDA B 

#$08 

431B 

CE 

00 

00 


LDX 

#$00 

431E 

FF 

42 

59 


STX 

H 

4321 

0C 



BRLP3 

CLC 


4322 

79 

42 

5A 


ROL 

L 

4325 

79 

42 

59 


ROL 

H 

4328 

46 




ROR A 


4329 

24- 

07 



BCC' 

BRLP4 

432B 

FE 

42 

59 


LDX 

H 

432E 

08 




I NX 


432F 

FF 

42 

59 


STX 

H 

4332 

5A 



BRLP4 

DEC B 


4333 

26 

EC 



BNE 

BRLP3 

4335 

F6 

42 

52 


LDA B 

LL8 

4338 

B6 

42 

5B 


LBA A 

D 

433B 

0C 



BRLP5 

CLC 


433C 

79 

42 

5A 


ROL 

L 

433F 

79 

42 

59 


ROL 

H 

4 342 

46 




ROR A 


4343 

24 

07 



BCC 

BRLP6 

4345 

FE 

42 

59 


LDX 

H 

4348 

08 




INX 


4349 

FF 

42 

59 


STX 

H 

434C 

5A 



BRLP6 

DEC B 


434D 

26 

EC 



BNE 

BRLP5 

434F 

39 




RTS 


4350 

01 



LLLE8 

NOP 


4351 

CE 

42 

0E 


LDX 

#LL 

4354 

B6 

42 

5C 


LDA A 

E ■ . 

4357 

E6 

00 



LDA B 

aX 

4359 

CE 

00 

00 


LDX 

#$00 

43 5 C 

FF 

42 

59 


STX 

H 

435F 

0C 



BRLPl 

CLC 


4360 

79 

42 

5A 


ROL 

■ L, 

4363 

79 

42 

59 


ROL 

H 

4366 

46 




ROR A 


4367 

24 

0 7 



BCC 

BRLP2 

4369 

FE 

42 

59 


LDX 

H 

436C 

08 




INX 


436D 

FF 

42 

59 


STX 

H 

4370 

5A 



BRLP2 

DEC B 


4371 

26 

EC 



BNE 

BRLPl 

4373 

39 





iliil® 


♦ MAIM, FFT PROGRAM STARTS 

♦ N HAS #QF SAMPLE POINTS 



4374 

4377 

437D 

4380 

4383 

4386 

4389 

438C 

438F 

4390 
4393 
4396 
4399 
4390 
439F 
43A2 
43A4 

43A5 
43A7 
4 3A8 
43AB 
43AD 
43AF 
43B2 
43B3 
43B5 
43B8 
43B9 
43BC 
43BE 
43BF 
43C0 
43C2 
43C3 
43C4 

■tm$ 

43C7 


43CA 
43 CD 
43D0 
43D3 
43D6 
43D9 
43DC 
43DF 
43E2 
4iSl 
43E8 
43EB 
43EE 


FE 

42 

0A 

FFT 

LDX 

N 

7F 

42 

59 


CLR 

H 

If 

m 

SH 



L 

7F 

42 

5B 


CLR 

D 

7F 

42 

5C 


CLR 

E 

7F 

42 

5E 


CLR 

C 

7F 

42 

5D 


CLR 

Bl- 

FF 

42 

59 


STX 

H 




* 

TO CALCULATE L0G2N 

BD 

51 

CC 


JSR 

L0G2N S 

0C 




CLC 


76 

42 

59 


ROR 

H 

76 

42 

5A 


ROR 

L 

FE 

42 

59 


LDX 

H 

FF 

42 

0C 


STX 

lA 

FF 

‘•■rfc L.! 

06 


STX 

NI2 

F6 

42 

5E 


LDA B 

C 

86 

08 



L DA A 

#S08 

11 




CBA 





* 

CHECKING L0G2N>8 OR <8 

2B 

06 



BMI 

GTS 

4F 




CLR A 


B7 

42 

51 


ST A A 

BFLAG 

20 

0B 



BRA 

LT8 ■ 

86 

01 


GTS 

LDA A 

#$01 

B7 

42 

51 


STA A 

BFLAG 

17 




TBA 


80 

08 



SUB A 

#$08 

B7 

42 

52 


STA A 

LL8 

17 



LT8 

TBA 


B7 

42 

0E 


STA A 

LL 

80 

02 



SUB A 

#$02 

16 




TAB 


0C 




CLC 


86 

01 



LDA A 

#$01 

46 




ROR A 


46 



LOOP 

RO R A 


5 A 




IlSrC B 





... : B. |E 

LOOP 

B7 

42 

10 


STA A 

HH 




* 

CALCULATION OF DIFFERENT 




, ^ 

OFFSETS TO GENERATE 





SINE AND CO S I N E VALUES 

CE 

40 

04 


LDX 

#SIN J 

FF 

42 

: 5B 


STX 

::-'D ■ ^ 

B6 

42 

16 


LDA A 

C1024 ; 

F6 

42 

: 1 7 


LDA B 

Cl 024+1 

FB 

42 

5C 


ADD B 

E 

B9 

42 

5B 


ADC A 

D 

B7 

42 

; 14 


STA A 

S1024 

F7 

42 

; 15 


STA B 

SI 024+1 

CE 

40 

04 


LDX 

#SIN 

FF 

42 

1 5B 




B6 

42 

1 ID 


■ LDA A 

'XXii-iiill 

F6 

42 1C 


.LDA B 

jxfiilixXx!; 

BB 

42. 

5C 


ADD A 

E 



43F1 

F9 

42 

5B 


ABC B 

D 

43F4 

B7 

42 

19 


STA A 

SI512+1 

43F7 

F7 

42 

18 


STA B 

SI512 

43FA 

■ CE 

40 

04 


LDX 

#SIN 

43FD 

FF 

42 

5B 


STX 

D 

4400 

B6 

42 

13 


LDA A 

'CT512+1 

4403 

F6 

42 

12 


LDA B 

CT512 

4406 

BB 

42 

5C 


ADD A 

E • 

4409 

F9 

42 

5B 


ADC B 

D 

440C 

B7 

42 

IB 


STA A 

SN512+1 

440F 

F7 

• 42 

lA 


STA B 

SN512 

4412 

4F 




CLR A 


4413 

B7 

42 

08 


STA A 

I 

4416 

CE 

00 

01 


LDX 

#$0001 

4419 

FF 

42 

1 E 


STX 

IB 





♦ 

OUTER LOOP 





* 

PARAMETERS SETTING 

44 1C 

CE 

00 

00 

OUTER 

LDX 

#$0000 ; 

44 IF 

FF 

42 

20 


STX 

IC 

4422 

FF 

42 

26 


STX 

K 

4425 

FE 

42 

0C 


LDX 

lA 

4428 

FF 

42 

59 


STX 

H 

442B 

FF 

42 

22 


STX 

I D 

442E 

0C 




CLC 


442F 

79 

42 

5A 


ROL 

L 

4432 

79 

42 

59 


ROL 

H 

44 35 

FE 

42 

59 


LDX 

H 

4438 

FF 

42 

24 


STX 

I2D 





* 

INNER 

LOOP PARAMETERS 





* 

SETTING 


443B 

FE 

42 

0C 

INNER 

LDX 

I A 

443E 

FF 

42 

5B 


STX 

D 

4441 

FE 

•42 

20 


LDX 

IC 

4444 

FF 

42 

59 


STX 

H 





♦ 

POWER OF W CALCULATED 

4447 

01 



ICIA 

NOP 


4448 

0C 




CLC 


4449 

76 

42 

5B 


ROR 

D 

444G 

76 

42 

5C 


ROR 

E . . 

444F 

25 

08 



BCS 

OICIA 

4451 

76 

42 

59 


ROR 

H 

4454 

76 

42 

5A 


ROR 

;L 

4457 

20 

EE 



BRA 

ICIA 

4459 

0 1 



01 CIA 

NOP 


445A 

BD 

43 

08 


JSR 

BI TR 

445D 

4F 




CLR A 


44 5 E 

B6 

42 

10 


LDA A 

HH 

4461 

46 



I CLP 

ROR A 


4462 

81 

00 



CMP A 

#$00 

4464 

27 

08 



BEQ 

FINIS 

4466 

78 

42 

5A 


ASL 

L 

4469 

79 

42 

59 


ROL 

H 

446C 

20 

F 3 



BRA 

ICLP 

446E 

01 



FINIS 

NOP 



r liv 1 D wur 

%! If' ' - V' :-i#t « ^ 



446F 

01 


GL256 

NOP 


4470 

B6 

42 59 


LDA A 

H 

4473 

81 

01 


CMP A 

#$01 

4475 

25 

3A 


BCS 

NLT257 

4477 

81 

01 


CMP A 

#$01 

4479 

26 

07 


BNE 

GT256 

447B 

F6 

42 5A 

ACHEC 

LDA B 

L 

447E 

Cl 

01 


CMP B 

#$01 

4480 

25 

2F 


BCS 

NLT257 




* 

VALUES OF SINE AND 



. 

* 

FOB W>256 

4482 

01 


GT256 

NOP 


4483 

78 

42 5A 


ASL 

L 

4486 

79 

42 59 


ROL 

H 

4489 

FE 

42 59 


LDX 

H 

448C 

FF 

42 5B 


STX 

D 

448F 

4F 



CLR A 


4490 

B6 

42 15 


LDA A 

S1024+1 

4493 

F6 

42 14 


LDA B 

S1024 

4496 

B0 

42 5C 


SUB A 

E 

4499 

F2 

42 5B 


SBC B 

D 

449C 

B7 

42 5A 


STA A 

L . 

449F 

F7 

42 59 


STA B 

H 

44A2 

FE 

42 59 


LDX 

H 

44A5 

A6 

00 


LDA A 

^X 

44A7 

08 



INX 


44A8 

E6 

00 


LDA B 

vX 

44AA 

CE 

42 28 


LDX 

#SINE 

44AD 

A7 

00 


STA A 

>X 

44AF 

20 

02 


BRA 

NEXT 

44B 1 

20 

2F 

NLT257 

BRA 

LT25 7 

44B3 

08 


NEXT 

INX 


44B4 

E7 

00 


STA B 

\»X 

44B6 

B6 

42 19 


LDA A 

SI 5 12 + 1 

44B9 

F6 

42 18 


LDA B 

SI 5 12 

44 BC 

BB 

42 5C 


ADD A 

E 

44 BF 

F9 

42 5B 


ADC B 

D 

44C2 

B7 

42 5A 


STA A 

L 

44C5 

F7 

42 59 


STA B 

H 

44C8 

4F 



CLR A 


44C9 

FE 

42 59 


LDX 

H 

44 CC 

A6 

00 


LDA A 

.^X 

44CE 

08 



INX 


44CF 

E6 

0 0 


LDA B 

/ j.'x- ■ 

44 Dl 

CE 

42 2B 


LDX ' 

#C0SIN+1 

44 D4 

43 



COM A 


44D5 

53 



COM B 


44D6 

CB 

01 


add S' 

#$01 

44D8 

89 

00 


ADC A 

#$'00 

44DA 

E7 

00 


STA B 

>x 

44DC 

09 



DEX 


44DD 

A7 

00 


STA A 

J-X 

44 DF 

7E 

45 40 


JMP 

INI 




* 

VALUES 

OF SINE AN 





FOR W<=1 



44E2 01 LT257 NOP 



44E3 

78 

42 5A 

ASL 

L 

44E6 

79 

42 59 

ROL 

H 

44E9 

FE 

42 59 

LDX 

H 

44EC 

FF 

42 5B 

SIX 

D 

44EF 

CE 

40 04 

LDX 

i^-SIM 

44F2 

FF 

42 59 

SIX 

H 

44F5 

B6 

42 5A 

LDA A 

L 

44F8 

F6 

42 59 

LDA B 

H 

44FB 

BB 

42 5C 

ADD A 

E 

44FE 

F9 

42 5B 

ADC B 

D 

4501 

B7 

42 5A 

STA A 

L 

4504 

F7 

42 59 

STA B 

H 

4507 

FE 

42 59 

LDX 

H 

450A 

A6 

00 

L Da a 

..X 

450C 

08 


INX 


450D 

E6 

00 

LDA B 

..X 

450F 

FF 

42 59 

SIX 

H 

4512 

CE 

42 28 

LDX 

, #SIME 

4515 

A7 

00 

STA A 

^X 

45 17 

08 


INX 


4518 

E7 

00 

STA B 

jK 

451A 

4F 


CLR A 


451B 

B6 

42 IB 

LDA A 

SN512.+ 1 

451 E 

F6 

.42 lA 

LDA B 

SN5i2 

4521 

B0 

42 5C 

SUB A 

E 

4524 

F2 

42 5B 

SBC B 

D 

4527 

B7 

42 5A 

STA A 

L 

452A 

F7 

42 59 

STA B 

H 

45 2D 

FE 

42 59 

LDX 

H 

453l?l 

A6 

00 

LDA A 

>X 

4532 

08 


INX 


4533 

E6 

00 

LDA B 

^X 

4535 

FF 

42 59 

STX 

H 

4538 

CE 

42 2A 

LDX 

#COSIN 

453B 

A 7 

00 

STA A 

rX 

453D 

08 


INX 


453E 

E7 

.00 

STA B 


4540 

01 

INI 

NOP 


4541 

FE 

42 0C 

LDX 

I A 

4544 

FF 

42 59 

STX 

H 

4547 

78 

42 5A 

ASL 

L : 

454A 

79 

42 59 

ROL 

H 

45 4 D 

FE 

42 59 

LDX 

: H 

4550 

FF 

42 5B 

STX 

D 

4553 

CE 

49 CA 

LDX 

#XREAL 

4556 

FF 

42 59 

STX 

. H 

4559 

B6 

42 5 A 

LDA A 

L' . 

455C 

F6 

42 59 

LDA B 

H 

455F 

BB 

42 5C 

ADD A 

E 

4562 

F9 

42 5B 

ADC B 

D 

4565 

B7 

42 39 . 

STA A 

MIE+1 

4568 

F7 

42,38 

STA B 

MIR 

456B 

CE 

4D CA. 

LDX 

#XIMA@ - 

456E 

FF 

42 59 


H 

4571 

B6 


' LDA A 


4574 

F6 

42 59 

. LDA B 




4577 

BB 

42 

5C 

ADD A 

E 

457A 

F9 

42 

5B 

ADC B 

D 

457D 

B7 

42 

3B 

STA A 

MII + 1 

4580 

F7 

42 

3A 

STA B 

Mil 

4583 

FE 

42 

20 

LDX 

IC 

4586 

FF 

42 

59 

SIX 

H 

4589 

78 

42 

5A 

ASL 

L 

458C 

79 

42 

59 

ROL 

H 

458F 

B6 

42 

5A 

LDA A 

L 

4592 

F6 

42 

59 

LDA B 

H 

4595 

B7 

42 

5E 

STA A 

C 

4598 

F7 

42 

5D 

STA B 

B1 




* 

INNER MOST LOO 

459B 

01 


IMOST 

NOP 


459C 

B6 

42 

5E 

LDA A 

C 

459F 

F6 

42 

5D 

LDA B 

B 1 

45A2 

36 



PSH A 


45A3 

37 



PSH B 


45A4 

BB 

42 

39 

ADD A 

MIR+1 

45A7 

F9 

42 

38 

ADC B 

MIR 

45AA 

B7 

42 

5A 

STA A 

L 

45AD 

F7 

42 

59 

STA B 

H 

45B0 

FE 

42 

59 

LDX 

H 

45B3 

A6 

00 


LDA A 


45B5 

08 



INX 


45B6 

E6 

00 


LDA B 

^X 

45B8 

FF 

42 

5B 

STX 

D 


B7 

42 

59 

STA A 

H 

m 45BE 

F7 

42 

5A 

STA B 

L 

45C1 

BD 

42 

E0 

JSR 

OFL 

45C4 

FE 

42 

59 

LDX 

H 

45C7 

FF 

42 

40 

STX 

XYRMI 

45CA 

FF 

42 

5B 

STX 

D 

45 CD 

FE 

42 

2A 

LDX 

COSIN 

45D0 

FF 

42 

59 

STX 

H 

45 D3 

BD 

42 

5F 

JSR 

DPMUL 

45D6 

FE 

42 

59 

LDX 

H 

45D9 

FF 

42 

48 

STX 

REAL 1 

^ 45DC 

33 



PUL B 


45DD 

32 



PUL A 


45DE 

BB 

42 

3B 

ADD A 

Mll+i 

45E1 

F9 

42 

3A : 

ADC B 

Mil 

45E4 

B7 

42 

5A 

STA A 

: 'L: 

45E7 

F7 

42 

59 . 

STA B 

H 

45EA 

FE 

42 

59 

LDX 

H 

45ED 

A6 

00 


LDA A 

; v X; 

45EF 

08 



INX 


45F0 

E6 

00 


LDA B 

>X 

45F2 

FF 

42 

5B 

STX 

D 

45F5 

B7 

42 

59 

STA A 

H 

45F8 

F7 

42 

5 A ■ : 

STA B 

L 

45 FB-'' 

BD 

42 

E0 

JSR 

OFL 

4ife: ^ 

FE 

42 

59 

LDX 


mm 

FF 

42 

42 

Stx 

XYIMI 

m 04 

FF 

42 

5B 


: 

4607 

f£ 

42 

28 

■ LDX 

SINE 

l^fi' 






460A 

FF 

42 59 

STX 


H 

460D 

B6 

42 5E 

LDA 

A. 

C 

4610 

F6 

42 5D 

LDA 

B 

B1 

4613 

36 


'PSH 

A 


4614 

37 


PSH 

B 


4615 

BD 

42 5F 

JSR 


DPMUL 

46 18 

FE 

42 59 

LDX 


H 

46 IB 

FF 

42 5B 

STX 


D 

46 IE 

B6 

42 49 

LDA A 


REAL 1 + 1 

46 21 

F6 

42 48 

LDA 

B 

REAL 1 

4624 

BB 

42 5C 

ADD A 


E 

4627 

F9 

42 5B 

ADC B 


D 

462A 

B7 

42 4B 

STA 

A 

BREAL+1 

462D 

F7 

42 4A 

STA 

B 

BREAL 

46 30 

FE 

42 42 

LDX 


XYIMI 

4633 

FF 

42 5B 

STX 


D 

4636 

FE 

42 2A 

LDX 


COSIN 

4639 

FF 

42 59 

STX 


H 

463C 

BD 

42 5F 

JSR 


DPMUL 

46 3F 

FE 

42 59 

LDX 


H 

4642 

FF 

42 4C 

STX 


I M AG 1 

4645 

FE 

42 40 

LDX 


XYRMl 

4648 

F F 

' 42 5B 

STX 


D 

46 4B 

FE 

42 -28 

LDX 


SINE 

464E 

FF 

42 59 

STX 


H 

4651 

BD 

42 5F 

JSR 


DPMUL 

46 54 

FE 

42 59 

LDX 


H 

46 57 

BD 

42 CB 

JSR 


COMPL 

46 5A 

B6 

42 5A 

LDA 

A 

L , 

465D 

F6 

42 59 

LDA 

B 

H 

4660 

BB 

42 4D 

ADD A 


IMAGT+1 

4663 

F9 

42 4C 

ADC B 


IMAGl 

4666 

B7 

42 5A ■ 

STA 

A 

L 

4669 

F7 

42 59 

STA 

B 

H 

46 6 C 

FE 

42 59 

LDX 


H 

46 6F 

FF 

42 4E 

STX 


BIMAG 

4672 

33 


PUL B 



46 73 

32 


PUL A 


4674 

B7 

42 5E 

STA A 

C 

46 77 

F7 

42 5D 

STA 

B 

B1 

46 7A 

CE 

49 CA 

LDX 


#XREAL 

467D 

FF 

42 59 

STX 


H 

4680 

BB 

42 5 A 

ADDA 

'L;^; . 

468 3 

F9 

42 59 

ADC B 

■H 

4686 

B7 

42 5A 

STA 

A 

L - 

4689 

F7 

42 59 

STA 

B 


468C 

FE 

42 59 

LDX 



468F 

A 6 

0 0 

LDA 

A 

jX 

4691 

08 


I NX 



4692 

E6 

00 

LDA 

B 

/X 

4694 

B7 

42 59 

STA A 


H 

4697 

F7 

42 5A 

STA 

B 

L 

46 9A 

BD 

42 E0 

JSR 


OFL 

469D 

FE 


LDX 


H 

46A0 

FF 

42 44 



AREAL 


CE 

4D CA 

LDX 


#XTMAG 



46A6 

FF 

42 59 

STX 

H 

46A9 

B6 

42 5E 

LDA A 

C 

46AC 

F6 

42 5D 

LDA B 

B1 

46AF 

BB 

42 5A 

ADD A 

L 

46B2 

F9 

42 59 

ADC B 

H 

46B5 

B7 

42 5A 

STA A 

L 

46B8 

F7 

42 59 

STA B 

H 

46BB 

FE 

42 59 

LDX 

H 

46BE 

A 6 

00 

LDA A 

vX 

46C0 

08 


INX 


46 Cl 

E6 

00 

LDA B 

^X 

46C3 

B7 

42 59 

STA A 

H 

46C6 

F7 

42 5A 

STA B 

L 

46C9 

BD 

42 E0 

JSR 

OFL 

46CC 

FE 

42 59 

LDX 

H 

46CF 

FF 

42 46 

STX 

AIMAG 

46D2 

B6 

42 4B 

LDA A 

BREAL+1 

46 D5 

F6 

42 4A 

LDA B 

BREAL 

46 D8 

BB 

42 45 

ADD A 

AREAL +1 

46 DB 

F9 

42 44 

ADC B 

AREAL 

46 DE 

B7 

42 5C 

STA A 

E 

46 El 

F7 

42 5B 

STA B 

D 

46E4 

CE 

49 CA 

LDX 

#XREAL 

46E7 

FF 

42 59 

STX 

H 

46 EA 

B6 

42 5A 

LDA A 

L 

46ED 

F6 

42 59 

LDA B 

H 

46F0 

BB 

42 5E 

ADD A 

C 

46F3 

F9 

42 5D 

ADC B 

BI 

46F6 

B7 

42 5A 

STA A 

, L ■ ■ ' 

46 F9 

F7 

42 59 

- STA B 

H 

46FC 

FE 

42 59 

LDX 

H 

46 FF 

B6 

42 5B 

LDA A 

D. . 

4702 

F6 

42 5C 

LDA B 

E- 

47 05' 

A7 

00 

STA A 

-rX 

4707 

08 


INX 


4708 

E7 

00 

STA B 

>x 



♦ 

A REAL 

-B REAL CALCULATED 

470A 

FE 

42 4A 

LDX 

BREAL 

470D 

FF 

42 59 

STX 

H 

4710 

BD 

42 CB 

JSR ■ 

COMPL 

47 1 3 

B6 

42 5A 

LDA A 

'L 

4716 

F6 

42 '59 

LDA B 

H:"'': ■■ 

4719 

BB 

42 45 

ADD A 

AREAL+I 

47 IG 

F9 

42 44 

ADC B 

AREAL 

471F 

B7 

42 5C 

STA A 

E ■ ■ 

4722 

F7 

42 5B 

STA B 

■ ■ D- ' 

4725 

B6 

42 5E 

LDA A 

c:-. ■ 

4728 

F6 

42 5D 

LDA B 

'■;Bi ■; 

472B 

BB 

42 39 

ADD A 

MIR+l 

472E 

F9 

42 38 

ADC B . 

, MIR 

4731 

B7 

42 5A 

STA A 

L 

4734 

F7 

42 59 

STA B 

H . 

4737 

FE 

42 59 

LDX 

H 

473A 

B6 

42 5B 


D ■ 

473D 

A7 

00 

STA A 

^x 

473F 

08 


INX 




4740 

F6 

42 50 

LDA B 

E 

4743 

E7 

00 

STA B 

j,X 



* 

AIMAG+BIMAG OALCULATED 

4745 

B6 

42 4F 

LDA A 

BIMAG+1 

4748 

F6 

42 4E 

LDA B 

BIMAG 

474B 

BB 

42 47 

ADD A 

AIMAG+1 

474E 

F9 

42 46 

ADO B 

AIMAG 

4751 

B7 

42 50 

STA A 

E 

4754 

F7 

42 5B 

STA B 

D 

47 57 

OE 

4D OA 

LDX 

#XIMAG 

47 5A 

FF 

42 59 

STX 

H 

475D 

B6 

42 5A 

LDA A 

L 

4760 

F6 

42 59 

LDA B 

H 

4763 

BB 

42 5E 

ADD A 

0 

4766 

F9 

42 5D 

ADO B 

B1 

4769 

B7 

42 5A 

STA A 

L 

476C 

F7 

42 59 

STA B 

H 

476F 

FE 

42 59 

LDX 

H 

47 72 

B6 

42 5B 

LDA A 

D 

4775 

A7 

00 

STA A 

j X 

4777 

08 


INK 


47 78 

F6 

42 50 

LDA B 

E 

47 7B 

E7 

00 

STA B 

^ X * 

477D 

FE 

42 4E 

LDX 

BIMAG 

4780 

FF 

42 59 

STX 

H 

478 3 

BD 

42 OB 

JSR 

COMPL 

4786 

FE 

42 59 

LDX 

H 

4789 

FF 

42 5B 

STX 

D 

478C 

B6 

42 4 7 

LDA A 

AIMAG+1 

478F 

F6 

42 46 

LDA B 

AIMAG 

4792 

BB 

42 50 

ADD A 

E , , 

4795 

F9 

42 5B 

ADC B 

D 

4798 

B7 

42 50 

STA A 

E 

479B 

F7 

42 5B 

' STA B 

D 

479E 

B6 

42 3B 

LDA A 

MIT + 1 

47 A 1 

F6 

42 3A 

LDA B 

Mil 

47A4 

BB 

42 5E 

ADD A 

0 

47A7 

F9 

42 5D 

ADC B 

B1 

47AA 

'B7 

42 5A 

■ STA A 

L 

47AD 

F7 

42 59 

STA B 

H . ■' 

47B0 

FE 

42 59 

LDX 

H 

47B3 

B6 

42 SB 

LDA A 

■ D 

47B6 

A 7 

00 

STA A 


47B8 

08 


INX 


47B9 

F6 

42 50 

LDA B 

E ■ 

47BC 

E7 

00 

STAB 


47BE 

FE 

42 5 D 

LDX 

:B1 ; : 

4 701 

08 


INX 


4702 

08 


INX 


4703 

FF 

42 5D 

STX 

B1 

4706 

B6 

42 5E 

LDA A 

0 

4709 

OE 

42 25 

LDX 

#I2D+1 

4700 

A1 

00 

CNP A 


47 OE 

^28; 

02 


IMOSTl 

47D0 

20 

03 

BIIA 

NEXT3 


NEXT' ITERATION OF 







* 

INNER MOST LOOP 

47D2 

7E 

45 

9B 

IMOSTl 

JMP 

IMOST 

47D5 

09 



NEXT 3 

DEX 


47D6 

B6 

42 

5D 


LDA A 

El 

47D9 

A1 

00 



CMP A 

jX 

47DB 

26 

02 



BNE 

IM0ST2 

47DD 

20 

03 



BRA 

NEXT4 

47DF 

7E 

45 

9B 

IM0ST2 

JMP 

IMOST 






NEXT ITERATION ^ 





* 

INNER LOOP 

47E2 

FE 

42 

0C 

NEXT 4 

LDX 

lA 

47E5 

FF 

42 

59 


STX 

H 

47E8 

78 

42 

5A 


ASL 

L 

47EB 

79 

42 

59 


ROL 

H 

47EE 

FE 

42 

59 


LDX 

H 

47F1 

FF 

42 

5B 


STX 

D 

47F4 

B6 

42 

21 


LDA A 

IC+1 

47F7 

F6 

42 

20 


LDA B 

IC 

4 7 FA 

BB 

42 

5C 


ADD A 

E 

47FD 

F9 

42 

5B 


ADC B 

D 

4800 

B7 

42 

21 


STA A 

IC+1 

48 03 

F7 

42 

20 


STA B 

IC 

4806 

B6 

42 

23 


LDA A 

ID+1 

4809 

F6 

42 

22 


LDA B 

ID 

480C 

BB 

42 

5C 


ADD A 

'E 

480F 

F9 

42 

5B 


ADC B 

D 

48 12 

B7 

42 

23 


STA A 

ID+1 

48 15 

F7 

42 

22 


STA B 

I D 

4818 

48 




ASL A 


4819 

59 




ROL B 


48 lA 

B7 

42 

25 


STA A 

I2D+1 

48 I D 

F7 

42 

24 


STA B 

I2D, 

4820 

FE 

42 

26 


LDX 

K 

4823 

08 




INX 


4824 

FF 

42 

26 


STX 

K 

4827 

FF 

42 

; 5B 


STX 

D 

482A 

CE 

42 

IF 


LDX 

#IB+1 

48 2D 

B6 

42 

5C 


LDA A 

E 

48 30 

A1 

00 



CMP A 

^X 

4832 

26 

02 



BNE 

INNERl 

4834 

20 

03 



BRA 

NEXT5 

48 36 

7E 

44 

1 3B 

INNER! 

JMP 

INNER 

4839 

09 



NEXT 5 

DEX 


48 3A 

B6 

42 

5B 


LDA A 

D 

483D 

A1 

00 



CMP A 

' ■•'X' 

483F 

26 

02 



BNE 

INNER2 

48 41 

20 

03 



BRA 

MEXT6 

4843 

7E 

44 

3B 

INNER2 

JMP 

INNER 





■ 

NEXT ITERATION 





^ , 

OUTER LOOP 

4846 

4F 



NEXT6 

CLR A 


4847 

B6 

42 

IF 


LDA A 

IB+1 

48 4A 

F6 

42 IE 


LDA B 

IB 

484D 

49 




ROL A 


484E 

59 




ROL B 


V 48 4F. 

B7 

42 

IF 


STA A 

IB+1 



4852 

F7 

42 

IE 

4855 

4F 



4856 

B6 

42 

0D 

48 59 

F6 

4 2 

0C 

485C 

56 



48 5D 

46 



485E 

B7 

42 

0D 

4861 

F7 

42 

0C 

4864 

7C 

42 

08 

4867 

B6 

42 

08 

48 6A 

CE 

42 

0E 

486D 

A1 

00 


486F 

26 

02 


4871 

20 

03 


48 73 

7E 

44 

1C OUTERl 

• 



♦ 

4876 

CE 

00 

00 NEXTE 

4879 

FF 

42 

5D 

487C 

B6 

42 

5E XLPl 

48 7F 

F6 

42 

5D 

4882 

B7 

42 

5A 

4885 

F7 

42 

59 

4888 

36 



4889 

37 



48 8 A 

BD 

4 3 

08 

488D 

B6 

42 

59 

4890 

B1 

42 

5D 

4893 

2B 

04 


4895 

27 

05 


4897 

20 

06 


4899 

7E 

49 

9E NLPCON 

489C 

7E 

49 

9 1 NCHECA 

489F 

01 


EXCHG 

48A0 

78 

42 

5A 

48A3 

79 

42 

59 

48A6 

FE 

42 

59 

48 A9 

FF 

42 

57 

48AC 

FF 

42 5B 

48AF 

78 

42 

5E 

48 B2 

79 

42 

5D 

48B5 

CE 

49 

CA 

48 B8 

FF 

42 59 

48BB 

B6 

42 

5A 

48BE 

F6 

42 

59 

48C1 

BB 

42 

5C 

48 C 4 

F9 

42 

5B 

48C7 

B7 

42 

5A 

48 C A 

FT 

42 

59 

48 CD 

FE 

42 

59 ■ 

48 D0 

A6 

00 


48 D2 

08 



48D3 

E6 

00 


48D5 

09 



48 D6 

FF 

42 5B 

48D9 

B7 

42 

54 

48DC 

F7 

42 53 


STA B IB 

CLR A 

LDA A IA+1 

LDA B I A 

ROR B 
ROR A 

STA A IA+1 

STA B I A 

INC I 

LDA A I 

LDX #LL 

CMP A 

BNE OUTER 1 

BRA NEXTE 

UMP OUTER 

UNSCRAMBLING OF OUTPUT 


LDX 


#$0000 

STX 


B1 

LDA 

A 

C 

LDA 

B 

B1 

STA 

A 

L 

STA 

B 

H 

PSH 

A 


PSH 

B 


JSR 


BITR 

LDA 

A 

H 

CMP 

A 

B1 

BMI 


NLPCON 

BEQ 


"NCHECA 

BRA 


EXCHG 

JMP 


LPCON 

JMP 


CHECA 

NOP 

ASL 


L 

ROL 


H 

LDX 


H 

STX 


BITRN 

STX 

D 

ASL 


C 

ROL 


B1 

LDX 


#XREAL 

STX 

H 

LDA 

A 

L 

LDA 

B 

■ H :: 

ADD 

A 

E 

ADC 

B 

: D 

STA 

A 

L 

STA 

B 

H 

LDX 


K 

LDA 

A 

,X 

INX 

LDA 

B 

^X 


DEX 

STX D 

STA A TREAL+1 
STA B tREAL 



48DF 
48E2 
48E5 
48E8 
48EB 
48 EE 
48FI 
48F4 
48F7 

48 FA 
48FG 
48FD 
48FF- 
4900 
4903 
4906 

4908 

4909 

49 0B 
490E 
49 1 1 
49 14 
4917 
4919 
49 lA 
491D 
491F 
4922 
4925 
4928 
49 2B 
492E 
49 3 1 
49 34 
49 37 
49 3A 
493D 
4940 

4942 

4943 

4945 

4946 
4949 
49 4C 
49 4F 
49 52 
4955 
49 58 
49 SB 
495E 
4961 
4964 
4967 
49 6A 
496C 
496D 


CE 

49 

CA 

LDX 

#XREAL 

FF 

42 

1 59 

SIX 

H 

B6 

42 

5A 

LDA A 

L 

F6 

42 

1 59 

LDA B 

H 

BB 

42 

5E 

ADD A 

C 

F9 

42 

5D 

ADC B 

B1 

B7 

42 

5A 

STA A 

L 

F7 

42 

59 . 

STA B 

H 

FE 

42 

59 

LDX 

H 

A6 

00 


LDA A 

^X 

08 



INX 


E6 

00 


LDA B 

^X 

09 



DEX 


FF 

51 

CA 

STX 

TEMP 

FE 

42 

5B 

LDX 

D 

A7 

00 


STA A 

,X 

08 



INX 


E7 

00 


STA B 

:.x 

FE 

42 

53 

LDX 

TREAL 

FF 

42 

59 

STX 

H 

FE 

51 

CA 

LDX 

TEMP 

B6 

42 

5A 

LDA A 

L 

A7 

00 


STA A 

/X 

08 



INX 


B6 

42 

59 

LDA A 

H 

A7 

00 


STA A 

^X 

FE 

42 

57 

LDX 

BITRN 

FF 

42 

5B 

STX 

D 

CE 

4D 

CA 

LDX 

#XIMAG 

FF 

42 

: 59 

STX 

H 

B6 

42 

5A 

LDA A 

L 

F6 

42 

59 

LDA B 

H 

BB 

42 

5C 

ADD A 

E 

F9 

42 

5B 

ADC B 

D 

B7 

42 

5A 

STA A 

L 

F7 

42 

59 

STA B 

H 

FE 

42 

59 

LDX 

H 

A6 

00 


LDA A 

^X 

08 



INX 


E6 

00 


LDA B 

>X 

09 



DEX 


FF 

42 

5B 

STX 

D 

B7 

42 

56 

STA A 

TIMAG + 1 

F7 

42 55 

STA B 

TTMAG 

CE 

4D 

CA 

LDX 

#XIMAG 

FF 

42 59 

STX 

: H : 

B6 

42 

5E : 

LDA A 

c 

F6 

42 5D 

LDA B 

B1 

BB 

42 

5A 

ADD A 

L 

F9 

42 

59 

ADC B 

H 

B7 

42 

5A 

STA A 

L 

F7 

42 

59 

STA B 

H 

FE 

42 

59 

LDX 

H 

A6 

00 


LDA A 


08 



INX 


E6 

00 


. LDA B 




496F 

09 




DEX 


4970 

FF 

51 

CA 


STX 

TEMP 

4973 

FE 

42 

5B 


LDX 

D 

4976 

A7 

00 



STA A 

^X 

4978 

08 




INX 


4979 

E7 

00 



STA B 


49 7B 

FE 

42 

55 


LDX 

TIMAG 

497E 

FF 

42 

59 


STX 

H 

4981 

FE 

51 

CA 


LDX 

TEMP 

4984 

B6 

42 

5A 


LDA A 

L 

4987 

A7 

00 



STA A 

#X 

4989 

B6 

42 

59 


LDA A 

H 

498C 

08 




INX 


498D 

A7 

00 



STA A 

>X 

498F 

20 

0D 



BRA 

LPCON 

499 1 

B6 

42 

5A 

CHECA 

LDA A 

L 

4994 

B1 

42 

5E 


CMP A 

C 

4997 

27 

05 



BEQ 

LPCON 

4999 

2B 

03 



BMI 

LPCON 

499B 

7E 

48 

9F 


JMP 

EXCHG 

499E 

' 33 



LPCON 

PUL E 


499F 

32 




PUL A 


49A0 

8B 

0 1 



ADD A 

#101 

49A2 

25 

02 



BCS 

NEXT8 

49A4 

20 

01 



BRA 

NEXT7 

49A6 

5C 



NEXT 8 

INC B 


49A7 

B7 

42 

5E 

NEXT 7 

STA A 

C 

49AA 

F7 

42 

5D 


STA B 

B1 

49AD 

FE 

42 

0A 


LDX 

N 

49B0 

FF 

42 

59 


STX 

H 

49B3 

B6 

42 

59 


LDA A 

H 

49B6 

B1 

42 

5B 


CMP A 

B1 

49B9 

27 

03 



BEQ 

SCHEC 

49 BB 

7E 

48 

7C 


JMP 

XLPl 

49BE 

■B6 

42 

5A 

SCHEC 

LDA A 

' L ■ 

49C1 

B1 

42 

5E 


CMP A 

C , 

49C4 

27 

0 3 



BEQ 

XEND 

49 C 6 

7E 

48 

7C 


JMP 

XLPl 





* 

END OF 

FFT SUBROUTINE 

49C9 

39 



XEND 

RTS 






* 

XREAL 

STORES 512 REAL VALUES 

49 CA 




XREAL 

RMB 

1024 ; 






XIMAG 

STO RES 5 12 I MAG • VALUE S 

4 DC A 




XIMAG 

RMB 

1024 ; 





* 

TEMP WORK LOCATIONS 

51CA 

00 



TEMP 

FCB 

00 L 

51 CB 

00 



TEMPI 

FCB 

/ 00: 





L0G2N 

L0G2N 

S U BRO UT I N E F IN D S LO G . OF 

5 ICC 

C6 

08 


LDA B 

#$08 

5 ICE 

B6 

42 

5A 


LDA A. 

L 

51D1 

0C 




CLC 


51D2 

46 



POWER 

ROR A 


51D3 

25 

0D 



BCS 

OUT 

51D5 

7C 

42 

SE, 


INC 

C . 



51D8 

5A 

' 

DEC B 


51D9 

26 

F7 

BNE 

POWER 

51DB 

C6 

08 

LDA B 

#$08 

51DD 

B6 

42 59 

LDA A 

H 

51E0 

20 

F 0 

BRA 

POWER 

51E2 

01 

OUT 

NOP 


51E3 

39 


RTS 


51E4 



END 



STATEMENTS =1184 
FREE BYTES =11986 


NO ERRORS DETECTED 



♦ ASSEMBLY nOl'IINt fOR 

♦ PSD MESlfREMFMT 

♦ ' 


49CA 




XREAL 

EQU 

i49CA 

ADCA 




XIMAG 

EQU 

S4DCA 

4BCA 




XREAL I 

ECU 

S4BCA 

AFCA 




XIMAG 1 

EQU 

i4FCA 

5600 




XRESUL^ 

EOU 

i5600 

425F 




DPMUL 

ECU 

$425F 

4259 




H 

EQU 

$4259 

425A 




L 

EQU 

$425A 

425B 




D 

EOU 

$425B 

425C 




E 

it 

EQU 

$425C 






SUBROUTINE OVERFLOW 

5400 





ORG 

$5400 

5400 

01 



OFL 

NOP 


5401 

78 

42 

5A 


ASL 

L 

5404 

79 

42 

59 


ROL 

H 

5407 

78 

42 

5A 


ASL 

L 

540A 

79 

42 

59 


ROL 

H 

540D 

78 

42* 

5A 


ASL 

L 

5410 

79 

42 

59 


ROL , 

H - 

5413 

78 

42 

5A 


ASL 

.. L^ 

5416 

79 

42 

59 


ROL 

H 

5419 

39 




RTS 






¥ 

♦ 

MAIN PSD ROUTINE STARTS 

5500 




¥ 

ORG 

$5500 

5500 

00 

00 


TEMP 

FDB 

$0000 

5502 

00 

00 


TEMPI 

FDB 

$0000 

5504 

00 

00 


TEMP2 

FDB 

$0000 

5506 

00 

00 


TEMP3 

FDB 

$0000 

5508 

00 

00 


TEMP4 

FDB 

$0000 

550A 

CE 

46 

CA 


LDX 

fXREALl 

550D 

FF 

55 

00 


STX 

TEMP 

5510 

CE 

4F 

CA 


LDX 

fXIMAGl 

5513 

FF 

55 

02 


STX 

TEMPI 

55 16 

CE 

56 

00 


LDX 

FXRESUL 

5519 

FF 

55 

08 


STX 

TEMP4 

55 1C 

FE 

55 

00 

INNER 

LDX 

TEMP 

55 IF 

A6 

00 


START 

LDA A 


5521 

08 




INK 


5522 

E6 

00 



LDA B 

' >X ; ■ 

5524 

08 




INX 


5525 

■FF- 

55 

, 00 


STX 

TEMP 

5528 

B7 

48 

59 


■:STa:::';a:'- 

H 

S52B 

F7 

4S 

r SA 


■■::5TA' B 

L 

SSBE'-" 

B0 

:'54 , 

00', 


JSR 

OFL 

5531 

■ FE 

42 

59 


LDX 

H 

5534 : 

rF,: 

42 

5t ^; ' 

A 

5Tx:; 

D 






DPMUL 

IS 16BIT MULTIPLY 



» jROUTINE IM FFT ROUTINF 


5537 

BB 

42 

5F 


aSR 

DPMUL 

55 3A 

EE 

42 

59 


L0« 

H 

5S3D 

rr 

55 

04 


STX 

TEMP2 

SS40' 

FE 

5S 

02 


LDX 

TEMPI 

5543 

A6 

00 



LDA A 

#X 

5545 

08 




INX 


SS46 

E6 

00 



LDA B 

#x 

S548 

08 




INX 


$549 

ff 

55 

02 


STX 

TEMPI 

SS4C 

B? 

42 

59 


STA A 

M 

5S4F 

¥7 

42 5A 


STA B 

L 

5552 

BO 

54 

00 


JSR 

OFL 

S555 

FE’ 

42 

59 


LDX 

H 

S5S8 

FF 

42 

SB 


STX 

D 

S55B 

3D 

42 

Sf 


aSR 

DPMUL 

SS5E 

FE 

42 

59 

, 

LDX 

H 

5561 

FF 

55 

06 


STX 

TEMPS 

5564 

Bd 

55 

05 


LDA A 

TEWPi^ 

5567 

r6 

55 

04 


LDA B 

TEMPB" * 

556A 

BB 

55 

07 


ADD A 

TEMP3-*-! 

556D 

F9 

55 

06 


ADC B 

TEMPS 

557® 

FE 

55 

08 


LDX 

TEMP4 

5573 

E7 

00 



STA B . 

^ #K. 

5575 

§6 




INK 


5576. 

A7 

00 



STA A 


S57S 

08 




INX 


5579 

FF 

55 

08 


STX 

TEMP 4 

557C 

FE 

55 

00 


LDX 

TEMP 

S57F 

8C 

4B 

GA 


CPX 

il<S4BCA 

5562 

87 

17 



: .. SEe 

EWDl 

5584 

8C 

40 

CA 


CPX 

#$4DCA 

5587 

27 

03 



BE® 

NEXT 

5589 

7E 

55 

JF 


- JMP 

START 

55 8C 

CE 

:49 

CA 

NEXT 

LDX 

iXREAL 

S$8F 

FF 

55 

0®' . 


STX 

TEMP 

5592 

CE 

40 

GA 


LDX 

#XiMAG 

5595 

FF 

55 

02 


STX 

TEMPI 

S598 

7E 

55 

1 C 


4MP 

INNER 

5598 

01 



ENDI 

: NOP 


5596 





: END 



STAtEHENtS ••98. : ; 
fBEE BYTES •24679 
m lEaORS BETECTED 



